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1. Introduction 

Collective phenomena are ubiquitous across physics and observed in systems having 
disparate temporal and spatial scales ranging from atomic nuclei, ultra-cold atomic 
gases and biological systems, up to the cosmological scales involved in galaxy for¬ 
mation. In a broad sense, the concept of collectivity refers to a common behavior 
exhibited by a group of entities (e.g., particles moving with a common velocity 
or along a preferred direction). Collectivity is often associated with an emergent 
phenomenon of a complex, many-body system, for which the basic interactions (for 
example, at the two-body level) may be well understood. An important goal of 
studying collective phenomena is to unravel how the macroscopic behavior of a 
many-body system emerges from its fundamental degrees of freedom. 

One of the most striking displays of collectivity is the strong expansion of the 
medium produced in ultra-relativistic heavy-ion collisions. A careful analysis of 
identified particle spectra, elliptic flow, Handbury-Brown-Twiss (HBT) interferom¬ 
etry and the suppression of particles at high transverse momenta has indicated 
that a hot and dense QCD medium, the “Quark-Gluon plasma (QGP)” undergoes 
a nearly ideal hydrodynamic expansion when created.^®^ 

Figure shows one of the many persuasive pieces of evidence for collective be¬ 
havior in heavy-ion collisions. The plot on the right shows the Arj-Ac/) correlation 
functions for pairs having 1 < px < 3 GeV/c in semi-central Pb-Pb collisions at 
the LHG. The right plot shows the same observable for minimum bias p-p colli- 
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Fig. 1. Two-particle correlation function for charged hadrons with transverse momenta 1 < pT < 
3GeV/c in minimum bias p-p collisions at ^/s = 7Te\^^ (left) and in 35-40% centrality Pb-Pb 
collisions at =2.76 Te\El (right), measured by the CMS experiment. 


sions where a dijet-like (or mini-jet) correlation structure can be seen; the narrow 
peak around (A 77 , A0) ^ (0,0) comes from the showering and hadronization of 
the leading parton, and the long-range away-side (A0 ^ tt) structure in A 77 rep¬ 
resents correlations from the recoiling parton. As the two partons are produced 
back-to-back they necessarily have opposite rapidities and will therefore populate 
the full acceptance in A 77 . In Pb-Pb collisions in addition to the jet-like correla¬ 
tions, a pronounced near-side (A(/) ^ 0 ) collimation extending over a long range 
in A 77 is observed and is now referred to as the “ridge”. This ridge-like correlation 
in heavy-ion collisions such as Au-Au at RHIC and Pb-Pb at the LHC is believed 
to be well understood: The overlap area of a heavy-ion collision at a finite impact 
parameter has an elliptic shape. The larger pressure gradients along the minor-axis 
of the ellipse lead to a larger flow in this direction and therefore collimated pro¬ 
duction in both directions of this axis creating a near- and away-side ridge. Careful 
subtraction of the away-side jet peak shows that such a double-ridge is present. The 
absence of a ridge structure in minimum bias p-p collisions suggest the absence of 
collective behavior in these systems even though event-by-event fluctuations may 
result in highly eccentric initial states. 

In 2010 surprising indications for collectivity in p-p collisions at the LHC were 
observed when triggering on rare events with high multiplicities (large number of 
final-state particles).The long-range near-side ridge for p-p events having charged 
multiplicity Atrk > 110 is shown in Fig. This ridge-like structure, not present 
in minimum bias p-p collisions, is reminiscent of the two-particle correlation in 
A-A collisions. Determining whether the p-p ridge can be attributed to collective 
flow effects will require a concerted effort by theorists and experimentalists. With a 
variety of theoretical proposals and limited experimental data in high-multiplicity 
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Fig. 2. Two-particle correlation function for particles having 1 < px < 3 GeV/c in high-multiplicity 
p-p collisions at ^/s = 7TeV, measured by the CMS experiment.!^ 


p-p collisions (see Ref.^^ for an early review) a conclusive explanation of the p-p 
ridge remains outstanding. 

With the discovery of the p-p ridge it was natural to look for a similar signal 
in p-A collisions. First collisions of high-energy proton and ion beams (p-A) were 
achieved at the LHC in 2011 at a center-of-mass energy of 5.02 TeV per nucleon 
pair. The discovery of a ridge in high-multiplicity p-Pb collisionJ^^tl^ did not come 
as a complete surprise, but the fact that the strength of the correlation was almost 
as large as the ridge in heavy-ion collisions was unexpected. 

The ridge correlations in p-p and p-A collisions remain to be fully understood, 
and the existence of similar structures in small colliding systems such as d- AiPH and 
^He-Au collision^SEll at lower RHIC energies have stimulated both experimental 
and theoretical communities to further investigate the properties of the ridge. 

This article provides a comprehensive review of the latest experimental re¬ 
sults and theoretical developments in our understanding of collective phenomena in 
dense, high-multiplicity hadronic systems. Perspectives on future directions are also 
discussed in light of future programs at major accelerator facilities such as RHIC 
and the LHC. 


2. Collectivity and hydrodynamics in A-A collisions 

Studies of multi-particle correlations have provided crucial insight into the underly¬ 
ing mechanism of particle production in relativistic heavy-ion collisions. The most 
prominent feature of multi-particle correlations in AA collisions is due to “elliptic 
flow”, an azimuthal anisotropy in momentum space induced by strong expansion 
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Fig. 3. Left: The 2-D and 1-D long-range Acf) two-particle correlation functions in ultra-central 
(0-0.2%) PbPb collisions at = 2.76 TeV measured by the CMS collaboration. Right: the 

extracted Vn Fourier harmonics from long-range two-particle correlations integrated over 0.3 < 
Pt < 3GeV/c, compared with hydrodynamic calculations.^^ 


of the initial almond-shaped overlap area of two nuclei.^ Elliptic flow phenom¬ 
ena can be well described by relativistic hydrodynamic simulations with a viscos¬ 
ity approaching the KSS bound r]fs = I/Att ^ 0.08.^^ In two-particle correlation 
measurements the elliptic flow generates an approximate cos(2A(/)) component to 
the two-particle correlation function that extends over a broad \Arj\ range .1^^ The 
long-range in rapidity is a consequence of the approximate boost invariance of the 
system. Studies of elliptic flow have been carried out over a wide range of energies 
and systemPHU at both RHIC and the LHC. 

Event-by-event fluctuations in the collision geometry of the initial-state intro¬ 
duce higher-order anisotropic flow components. Eor example, “triangular flow” re¬ 
sults in a cos(3A0) modulation of the two-particle correlation.^^^J^^ This contribu¬ 
tion can dominate the structure of the entire Acj) correlation in very central A-A 
collisions, as is shown in figure for the top 0-0.2% central Pb-Pb collisions at 
LHC. A prominent double hump appears on the away-side. 

By projecting the long-range component {\Ar]\ > 2) of the 2-D correlation func¬ 
tion onto a 1-D A0 axis the correlations can be characterized by a Eourier series, 
^1-1-2 v‘^ cos(nA0), where Vn denotes the single-particle anisotropy harmonic 

of order A dominant third-order Eourier component is seen because in ultra¬ 
central A-A collisions the impact parameter is small enough such that the initial 
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Fig. 4. Model calculations compared to measurements of the harmonic decomposition of azimuthal 
correlations produced in heavy ion collisions.^^The left panel shows model calculations and data 
for Vn vs. collision centrality in Pb+Pb collisions at sJsjsin =2.76 TeV. The right panel shows 
similar studies for the pT dependence of Vn in 200 GeV Au+Au collisions. The comparison of the 
two energies provides insight on the temperature dependence of 77 /s. 


(average) ellipticity is negligible. On an event-by-event basis, the initial geometry 
is governed by fluctuations resulting in comparable elliptic and triangular asymme¬ 
tries. 

Of course, any set of data can be decomposed into a Fourier series and the 
existence of a large elliptic anisotropy cannot be taken as proof of hydrodynamic 
behavior in and of itself. Quantitative comparison with theoretical calculations, in 
this case viscous relativistic hydrodynamic simulations is necessary in order to draw 
any strong conclusions. On the right-hand side of Fig. the extracted V 2 -Vq data 
in ultra-central Pb-Pb collisions from CMS are compared to viscous hydrodynamic 
calculations using two different initial condistions (the details of which will be 
discussed later). These results indicate that the produced medium has a shear 
viscosity to entropy density ratio ofry/s^O.OS — 0.2. 

Higher-order flow components in A-A collisions can provide more stringent con¬ 
straints on extracting both transport properties (e.g., r]/s) and initial-state models 
of heavy-ion collisions. Indeed, the elliptic and higher-order flow phenomena have 
been extensively studied at RHIC and the LHC over a wide range of collision 
centrality and particle px- Fig. shows the comparison of px-integrated Vn as a 
function of centrality in Pb-Pb collisions at the LHC (left) and Vn as a function of 
Px for 30-40% central An-An collisions at RHIC with hydrodynamic calculations 
using the IP-Glasma initial-state modelEHIMI The LHC data is well described by 
an 77 /s = 0.2 and at RHIC by 77 /s = 0.12 providing indication of a temperature de¬ 
pendence. It has now been widely accepted that in relativistic heavy-ion collisions, 
a strongly-coupled medium is formed exhibiting near-ideal fluid behavior. 
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3. Experimental evidence of collectivity in small systems 

This section provides a broad overview of experimental results in small colliding sys¬ 
tems, focusing on the recent results from p-Pb collisions at the LHC. While many of 
the experimental observables explored in this section will have striking similarities 
with the results from A-A collisions, one must be cautioned that strong conclu¬ 
sions about the interpretation of these measurements in A-A could only be made 
after detailed models (in this case hydrodynamics) confronted the data. Clearly, 
the system size in p-p or p-A collisions is significantly smaller, and therefore the 
applicability of hydrodynamics must be scrutinized. This section collects the wealth 
of experimental data we have on small systems and compares it to that in A-A with 
interpretations postponed to section]^ 

An important point that should be kept in mind throughout this discussion is 
that the nature of the very events we are studying in small systems are in a class of 
their own. They represent a fraction of the total cross-section, as can be gleamed 
from figure that shows the multiplicity distributions for minimum bias p-p, p-Pb 
and peripheral 50-100% central Pb-Pb collisions from CMS0IM1 Here, multiplicity 
is defined as the number of charged tracks within I 77 I < 2.4 and px > 0.4GeV/c. 
Events with 100-200 tracks (these are the high-multiplicity events where a ridge 
signal is seen in p-p and p-Pb) are a common occurrence in Pb-Pb. In general, 
experimental signatures should not depend only on the multiplicity; p-p and Pb-Pb 
events producing 100 tracks have a very different nature. 



Fig. 5. Multiplicity (A^trk) distributions for MinBias p-p, MinBias p-Pb and 50-100% centrality 
Pb-Pb collisions Ullill 
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3.1. Two-particle correlations and azimuthal anisotropy 

In September 2012, first collisions of proton and lead have been delivered by the 
LHC at = 5.02 TeV. About two million minimum bias p-Pb events were 

collected in a few hours of beam time by each experiment. Experiments with p-A 
were to serve as a control experiment for A-A collisions providing baselines for cold 
nuclear matter effect s.^^ 

However, with the recent discovery of the ridge in p-p, it was not clear what 
would be observed in high multiplicity p-A collisions. The p-A pilot run discovered 
a significant near-side ridge, as shown in Fig.[^ within the top few % (CMS,^ 
ATLA^ or even 20% (ALICEP) multiplicity class. A 2013 high-intensity p-Pb run 
provided access to higher multiplicity classes and increased statistics. 


pPb = 5.02 TeV at the LHC 



Fig. 6. The 2D two-particle correlation functions in high-multiplicity p-Pb collisions at = 

5.02 TeV measured by the ALICE (0-20%, left), El ATLAS (0-2%, middlej^ and CMS (0-3%, 

right jiSl 

experiments. 


The magnitude of the (near side) ridge in p-Pb collisions stands out as a promi¬ 
nent feature being about four times larger than the p-p ridge for >110 and 

becomes comparable to the away-side jet peak. In order to disentangle the ridge 
from the jet the ALICE and ATLAS collaborations subtracted the contribution of 
jets to the away-side correlation structure in high-multiplicity events, based on the 
estimates from low-multiplicity events. Implicit in this procedure is the assumption 
that there is negligible modification of the jet due to the high multiplicities (while 
this is clearly not valid in A-A collisions there has yet to be an observation of jet 
quenching in p-p or p-A). Following this subtraction procedure a “double” ridge 
structure having near- and away-side correlations of similar strength has been ob¬ 
served, as shown in Fig. for the 2D correlation function and its projection to ID 
in azimuthal angle. 

Motivated by the last decade of study of flow harmonics in A-A collisions the 
“double” ridge structure in p-Pb collisions has been analyzed using the same Fourier 
decomposition. The second- and third-order anisotropy Fourier harmonics, V 2 and 
Vs, are extracted from the long-range correlations as a function of px in high- 
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Fig. 7. The difference of two-particle correlation functions between high- and low-multiplicity p-Pb 
collisions at =5.02 TeV measured by the ALICE experiment.!^ 


multiplicity p-Pb collisions at = 5.02 TeV, shown in Fig. (bottom) from 

three LHC experimentsHHMEZI The curves shown from the CMS experiment are 
the results for V 2 and without subtracting the away-side jet demonstrating that 
in high-multiplicity events the contribution from jets to the Fourier harmonics are 
negligible at low enough px- 

Results of V 2 and in Pb-Pb collisions at = 2.76 TeV are also shown in 

Fig. 1^ (top), for a similar multiplicity range as for p-Pb collisions. ALICE, ATLAS 
and CMS results are in good agreement. Very recently, a sizeable V 2 has also been 
extracted from p-p collisions at 13 TeV by the ATLAS collaboration.!^ 

As shown in Fig. the V 2 and values first rise with px up to around 3 GeV 
and then fall off toward much higher px, a behavior that is very similar to Pb-Pb 
collisions. This may be an indication of a common origin of the ridge phenomenon 
in all collision systems. In Pb-Pb collisions, the trend at low px is consistent with 
hydrodynamic predictions, while the decrease of higher px Vn could be due to the 
lack of thermalization of more energetic probes. Surprisingly, the harmonics as 
a function of px for p-Pb and Pb-Pb systems are almost identical. This striking 
similarity extends to the full multiplicity range, as one can see from Fig. I (right) 
for V 2 , averaged over 0.3 < px < 3 GeV/c in p-Pb and Pb-Pb. Considering the 
drastically different initial-state geometry and system size of the two systems, this 
similarity and its implication remain to be fully understood. Arguments based on 
the universality of fluctuation-driven initial-state anisotropieJ^^®!! followed by the 
linear response of hydrodynamics may explain the coincidental triangular flows seen 
in p+Pb and Pb-Pb collisions. The V 2 value at a given multiplicity (Fig. left) is 
larger in Pb-Pb than in p-Pb; this is expected since we know in Pb-Pb there is a 
large initial eccentricity that after expansion generates a large V 2 . 
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Fig. 8. The second-order (y^) and third-order {y-^^ Fourier harmonics, extracted from long-range 
two-particle correlations, as a function of px" in Pb-Pb collisions at y'Xvw — TeV (top) and 
high-multiplicity p-Pb collisions at — 5.02 TeV (bottom). The curves represent original 

Vn data, while the markers denote the results after correcting for back- to-bac k jet correlations 
estimated from low-multiplicity events. Data are obtained from RefsHJllEZIl] 


3.2. Multi-particle azimuthal correlations 

One of the key questions about the nature of the ridge and its collectivity is whether 
it is only a two-particle correlation or if a correlation persists amongst all produced 
particles. The multi-particle cumulant technique was developed to measure the 
strength of collective azimuthal anisotropy by correlating four or more particles 
simultaneously .I33H311 Jt is designed to extract cumulants of collective azimuthal cor¬ 
relations, while suppressing non-collective short-range correlations, involving only 
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Fig. 9. The pT-averaged V 2 (left) and v^ (right) for 0.3<pt <3 GeV/c as a function of multiplicity 
in pPb (from two-particle correlation method) and PbPb (from event-plane method) collisions, 
measured by the CMS collaboration.!^ The solid line shows a hydrodynamic calculation using 
IP-Glasma initial conditions for Pb-Pb collisions.!^ 


a few particles. This approach has been widely used in studying A-A collisions. 


Measurements of the four-particle cumulant have been performed by the AT- 
LASp^ and ALICEP^ collaborations. It was later extended to six-, eight- 

and a//-particle (so-called “Lee-Yang Zeroes” (LYZ)) cumulants.l^ Fig. 10 shows 
the measured multi-particle cumulants, C 4 , ce, and cg, as a function of multiplicity 
in p-Pb and Pb-Pb collisions. These cumulants, Cn{n}, are calculated as follows: 


Cn{2}=((2)), 

Cn{4}=((4))-2.((2))^ 

Cn{6}=((6))-9.((4))((2)) + 12.((2))^ 

cn{8} =((8)) - 16 . ((6))((2)) - 18 . ((4))^ + 144 • ((4))((2))^ - 144((2))^ 


where two- and multi-particle azimuthal correlations are evaluated as: 

(( 2 )) = 

((6)) = ((e*’^('^l+'^2+03-04-05-06)^^^ 

= ^^g*’^(01+02+03+04-05-06-07-08)^^ 


( 2 ) 


Here ((•••)) represents the average over all combinations of particles from all events. 
The elliptic flow harmonics are related to the multi-particle cumulants by 

1'n{4} = yc„{4}, 1'n{6} = , Vn{8} = (3) 

which are shown in Fig. El for p-Pb and Pb-Pb systems side-by-side as a function 
of multiplicity. 
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Fig. 10. The four-(c 4 ), six-(c6) and eight-(c8) particle cumulants as a function of A^trk? averaged 
over (0.2)0.3 < pt < 3.0 GeV/c, in Pb-Pb collisions at = 2.76 TeV and p-Pb at = 

5.02 TeV. Data are obtained from RefslHMlI 


Within experimental uncertainties, V 2 values from various higher-order cumulant 
methods are all consistent with each other in both p-Pb and Pb-Pb as a function 
of multiplicity lending support to the highly collective nature of these systems. The 
V 2 from two-particle correlations does not follow the trends of the V 2 from higher 
cumulants. Part of the deviation may be due to the away-side jet which persists 
out to rapidity separations larger than the two unit rapidity gap imposed in the 
measurement. However, one would expect the effect from mini-jets to be smaller 
at higher multiplicities. The breaking of V2{‘2} ^ V 2 {^} is present in hydrodynamic 
models when event-by-event dynamical fluctuations of V 2 contribute differently to 
two-particle and higher-order cumulants. 

It should be stressed that although the observed scaling in Fig.[^is a necessary 
outcome of a hydrodynamic framework it is not sufficient proof that hydrodynamics 
is the correct underlying theory. The measurements are a milestone as they lend 
credence to the highly collective nature of p-Pb collisions and show that the ridge 
is a highly non-trivial emergent phenomena. 


3.3. Identified particle spectra and correlations 

Single-particle px spectra for various particle species as well as their multiplic¬ 
ity dependence in each collision system provide rich information for constraining 
the dynamics of particle production. Multiplicity dependence of identified particle 
(PID) px spectra have been measured in p-p, p-Pb and Pb-Pb collisions at the 
LHC.I^^t^In Fig. 12 (top), the average px {{pt)) for tt^/tt”, K+/K“, p/p. Kg and 
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Fig. 11. The second-order Fourier harmonic, r' 2 , as a function of A^trk obtained from two-, four-, 
six- and eight-particle cumulants, and the LYZ method, averaged over 0.3 < pt < 3.0 GeV/c, in 
Pb-Pb at = 2.76 TeV (left) and p-Pb at = 5.02 TeV (right).I^ 


A/A particles as measured by the ALICE collaboration is shown as a function of 
event multiplicity in p-Pb collisions at = 5.02 TeV.^ 

As multiplicity increases, px spectra of all particle species become flatter (or 
“harder”). A similar trend is also found in MC models such as PYTHIA and HI- 
JING, HIJING, which may be due to the biasing towards higher processes 
when triggering on high-multiplicity events. However, the most distinct feature of 
the data is not just the fact that (px) increases with multiplicity, but the increase 
is observed to be faster for particles with a heavier mass. For example, the ratio 
of proton to pion px spectra is enhanced in the higher px region, going from low- 
to high-multiplicity events. Moreover, the GMS collaboration reported that the 
amount of increase in (px) for each particle species from low to high multiplicity 
events is approximately proportional to the particle’s mass.l^ 

Similar mass dependence of PID spectra in A-A collisions has been extensively 
studied. In hydrodynamic models, it is attributed to the effect of a common velocity 
field along the radial direction of the medium, which generates a boost in particle 
momentum that is roughly proportional to the mass of the particle. A common 
framework for describing this characteristic mass dependence of px spectra is the 
“blast-wave” model,l^with parameters of a common kinetic freeze-out temperature, 
Tkin, and an average radial expansion velocity, (/3x) {Pt also depends on the radius). 
Fitting to the experimental spectra, (/3x) and Tkin can be extracted, and are 
(bottom) for p-Pb and Pb-Pb collisions for different multiplicity 


12 


shown in Fig. 

or centrality ranges. In Pb-Pb, more central collisions tend to have a larger (/dx) 
(due to larger pressure gradients) and a smaller T^in (because of a larger system size 
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Fig. 12. Top: Average transverse momentum of tt^/tt”, K+/K“, p/p, Kg and A/A particles as a 
function of multiplicity in pPb collisions at = 5.02 TeV. Bottom: parameters of blast-wave 

fits to identified particle spectra for pPb and PbPb collisions, and PYTHIAS models.I^ 


resulting in a longer lifetime of the hadronic rescattering stage). A similar trend is 
also observed in p-Pb collisions. Interestingly, if comparing p-Pb and Pb-Pb systems 
at similar multiplicities (about the same Tkin), the radial flow velocity is larger in 
the smaller p-Pb system. If fixing the total energy or entropy, a smaller QGP fluid 
possesses a stronger radial pressure gradient,I^^M^ which is consistent with the data. 
To verify this picture, one important step will be to carry out the same studies in 
high-multiplicity p-p collisions, which have an even smaller system size. 
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Fig. 13. The second-order Fourier harmonic, r’ 2 , for charged hadrons, tt^/tt”, K+/K“, and p/p 
from ALICE,^^ Kg and A/A from CMS, as a function of pT for 185 < A^trk < 220 in Pb-Pb 
collisions at = 2.76 TeV and p-Pb collisions at = 5.02 TeV. 


While the hydrodynamic fluid picture provides a natural interpretation to the 
data, as shown in Fig. (bottom), the PYTHIA model with color reconnec- 
tionJMUm can also generate similar behavior as seen in the data. In PYTHIA an 
outgoing hard part on and the underlying event are connected through color strings. 
As a result, strings will be pulled out of the underlying event in the transverse di¬ 
rection by the fast-moving parton, and then fragment into final-state particles. The 
color reconnection method increases this effect by reconnecting partons close in 
phase space. In this way, the “radially boosted” string provides an effective com¬ 
mon velocity for the hadrons produced from it, just like a fluid cell in the hydrody- 
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namic picture. Meanwhile, models based on gluon saturation can also qualitatively 
describe the general trend of the data.^^ 

To further differentiate theoretical models in p-Pb collisions, studies of identified 
particles have been extended to two-particle correlations. A long-range near-side 
two-particle correlation involving an identified particle is also observedEHlUS Results 
for the single-particle V 2 harmonic for j'K~ ^ K+/K“, p/p. Kg and A/A particles 
as a function of px in Pb-Pb and p-Pb collisions at similar multiplicities are shown 
in Fig. For both systems, a particle species dependence of is observed. At 
a given px, particles having a larger mass tend to exhibit a smaller anisotropy. 
This mass ordering effect of V 2 was first seen in A-A collisions at RHIC and LHC 
energiesplEllllHZ^ which can again be understood as the effect of radial flow pushing 
heavier particles toward higher-px pUESl Comparing p-Pb and Pb-Pb systems at 
similar multiplicities in Fig. a large mass splitting can be seen in the smaller p- 
Pb system, which is again consistent with the hydrodynamic pictureEHEIllt remains 
to be seen whether other theoretical interpretations, such as string fragmentation 
with color reconnections, can describe the features found in both single-particle and 
two-particle correlation data, as well as their system size dependence. 
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Fig. 14. Top: The second-order Fourier harmonic, for Kg, A/A, and unidentified charged 
particles as a function of px for 185 < Mrk < 220 obtained in p-Pb collisions at ~ 

5.02 TeV. Bottom: the rig-scaled values of Kg and A/A as a function of (KE^) jnq. Ratios of 
v^jnq to a smooth fit function of v^jnq for Kg as a function of (iFF^xf/rig are also shown.I^ 


In heavy ion collisions, a scaling phenomenon of V 2 with the number of con¬ 
stituent quarks (n^) has been discovered.^^^^^ The values of V 2 lnq are found to 
be very similar for all mesons {uq = 2) and baryons {uq = 3) when compared 
at the same transverse kinetic energy per constituent quark [KE^/riq^ where 
KE^ — V w? + Px — m). This empirical scaling may indicate that final-state 
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hadrons are formed through recombination of quarkspossibly providing evi¬ 
dence of deconfinement of quarks and gluons in these systems. This quark number 
scaling is found to be also valid within about 10% in high-multiplicity p-Pb colli¬ 
sions as shown in Fig. 14 for Kg and A/A particles. At similar multiplicities, the 


scaling holds better in the smaller p-Pb system than Pb-Pb system. Although the 
idea of quark coalescence is remarkably simple, the detailed dynamics of quark re¬ 
combination are far from fully understood. The observed constituent quark scaling 
in small systems may help elucidate this phenomena. 


3.4. Flow factorization 

When experimentally extracting the single-particle anisotropy Vn from multi¬ 
particle correlations the assumption is made that the multi-particle momentum 
distribution can be factorized into a product of single-particle distributions in each 
event, 

d^Nd^N d^N 

-=-• • •-. (4) 

dpidp2 • • • dpn dpi dp2 dpn 

The factorization assumption in equation 4 is broken by short-range correla¬ 
tions, for example, from jets, resonances, and momentum conservation. However, 
over a wide rapidity range (e.g., IA 77 I > 2 ) in A-A collisions, the predominant 
source of long-range correlations comes from an expanding fluid-like system with a 
common preferred flow direction (i.e., short-axis of the elliptical overlap region, or 
“event plane”). Other possible sources include back-to-back jet correlations on the 
away side but their contributions are negligible, especially for central A-A collisions, 
since they are strongly suppressed by the large multiplicity (1/Wrk) and also by 
the effect of jet quenching. Therefore, the relation in Eq.[^ holds and single-particle 
Vn can be extracted from the measurement of particle correlations. The factoriza¬ 
tion assumption has always been explicitly or implicitly applied in all methods of 
measuring anisotropic flow. Any breakdown of factorization would be an indication 
of correlations not originated from hydrodynamic flow. 

However, it was realized recently that even if all particles were to originate from 
a hydrodynamically flowing medium, the lumpy energy density distribution in the 
initial state would also break the factorization relation EEH This is because particles 
produced at different px and r] do not in general share a common flow direction, or 
event plane. The breakdown of the factorization relation as a function of particle’s 
Px and T] may provide information on the details of initial-state fluctuations, or 
more specifically, the lumpiness of the initial-state geometry. 

The px dependence of factorization breakdown has been quantified in terms of 
a factorization ratio, 

1%, A (Px ’ Pt ) ^nA (Px ’ Pt ) 



( 5 ) 
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Fig. 15. Elliptic flow factorization ratio, r 2 {p^,p^), as a function of - p^ , for 2.5 < p^ < 
3.0 GeV/c in ultra-central (0-0.2% centrality) Pb-Pb collisions at = 2.76 Te\^^ (left) and 

high-multiplicity p-Pb collisions at = 5.02 TeVEIEZIlll (right), compared to theoretical 

calculations in hydrodynamic modelslMMll 


where VnAiPT^Pr) represents the Fourier coefficient of two-particle azimuthal cor¬ 
relations for a given range. The data are shown in Fig. 15 for ultra-central 

(0-0.2% centrality) Pb-Pb (left) and high-multiplicity p-Pb (right) collisions. If the 
factorization relation holds, r 2 {p^,p^) will be constant at unity. Up to 20% devia¬ 
tion from unity of r 2 has been observed in ultra-central Pb-Pb collisions, which is 
consistent with hydrodynamic calculations using a Glauber initial condition, while 
the MC-KLN initial condition predicted a much larger effect .1^ This quantify is 
found to be insensitive to the r]/s value, providing very powerful constraints to 
models of initial-state fluctuations. Many follow-up theoretical studies showed that 
the magnitude of factorization breakdown is strongly correlated with the gran¬ 
ularity of initial-state fluctuationsElIimil A lumpy initial state with larger radial 
excitations would lead to a more significant breakdown of flow factorization, as seen 
in ultra-central Pb-Pb collisions. On the other hand, very little effect is found in 
high-multiplicity p-Pb collisions (Fig.[^ right). This may reflect a relatively larger 
granularity of initial-state fluctuations in p-Pb with respect to its smaller system 
size. Or in other words, a p-Pb system is relatively smoother than an ultra-central 
Pb-Pb system. Note that for both systems, there is almost no average geometry. 
Any final-state anisotropy is entirely generated by fluctuations. 


3.5. Femtoscopy from Hanbury-Brown-Twiss correlations 

Hanbury-Brown-Twiss (HBT) correlations provide crucial space time information 
of the source at the freezeout surface. The extracted femtoscopic radii in various 
collisions systems from p-p, p-A and A-A at RHIC and the LHC are summarized in 
Fig. 16 as a function of event multiplicity.!^ The femtoscopic radii follow a linear 
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1 /3 

^trk dependence for all systems, although the slopes are somewhat different. At 
similar multiplicities, the femtoscopic radii (extracted from a Gaussian fit to the 
two-particle correlation function) in p-Pb collisions is about 15-20% larger than that 
in p-p collisions. Both are significantly smaller than the value in Pb-Pb collisions, 
possibly due to a smaller initial system size. We will return to the discussion of 
HBT correlations when we present the details of hydrodynamic model calculations 
in later sections. 
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Fig. 16. Comparison of femtoscopic radii (extracted from a Gaussian fit to the two-particle 
correlation function), as a function of charged-particle multiplicity, measured for various collision 
systems and energies .1^ 


3.6. Heavy flavor production 

Production of heavy flavor quarks in p-p and p-Pb collisions has also been studied as 
a function of multiplicity. Studies of heavy quarkonia states (e.g., J/^, T) in heavy 
ion collisions can probe the possible onset of color screening from a deconfined QGP 
medium.l^ An intriguing multiplicity dependence of the T(2S) to T(1S) yield ratio 
from p-p, p-Pb to Pb-Pb collisions has been observed by the GMS collaboration.!^ 
The yield of the T(2S) state is suppressed relative to the T(1S) state as multiplicity 
increases in p-p and p-Pb collisions. Although limited in multiplicity reach, the 
amount of relative suppression has a tendency of smoothly extrapolating to the 
values measured in Pb-Pb collisions at higher multiplicities. This may be a hint 
that similar medium effects on the T states become significant in p-p and p-Pb 
collisions as multiplicity increases. However, more studies are still needed to exclude 
any possible bias introduced by the multiplicity selection. 

Open heavy flavor production has been measured in p-p collisions as a function 
of multiplicity by the ALIGE collaboration.!^ The yields of open charm 
and mesons as a function of multiplicity in p-p collisions at ^/s = 7TeV are 
shown in Fig. 17 Both the yield (y-axis) and multiplicity (x-axis) are normalized 


by the average values from minimum bias events. The open charm meson yields 
first linearly increase with multiplicity, which is consistent with the expectation of 
independent multiple parton interactions. At above 3-4 times the average minimum 
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bias multiplicity, the increase in yield becomes much faster than linear. Around a 
similar multiplicity region, the long-range ridge correlations also start becoming 
significant. These observations suggest the emergence of new dynamic processes in 
very high-multiplicity p-p collisions. 
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Fig. 17. Left: ratios of T(2S) to T(1S) yields as a function of charged-particle multiplicity in p-p 
collisions at ^/s = 7 TeV, p-Pb at = 5.02 TeV and Pb-Pb at = 2.76 TeV.I^ Right: 

average D meson and non-prompt J/T relative (to minimum bias events) yields as a function of 
the relative (to minimum bias events) charged-particle multiplicity at midrapidity in p-p collisions 
at = 7 TeV.II^ 


4. Theoretical interpretations 

Since the 2010 discovery of the ridge in p-p collisions by the CMS collaboration 
there has been a flurry of theoretical work to uncover its origin. Five years later, as 
the dust has started to settle, the heavy-ion community has ruled out and refined 
various models. In this section we will provide an overview of some of the theoretical 
ideas on the market, focusing on those that confront the data on a quantitative level. 

The dialogue in the heavy-ion community has been divided between initial-state 
and final-state models depending on whether the momentum space anisotropy is es¬ 
tablished at the moment of particle production or subsequently during strong final 
state interactions, or (while clearly not equivalent but sometimes used interchange¬ 
ably) non-hydro dynamic vs hydrodynamic descriptions. However, hydrodynamics 
is an effective description of the long wavelength, small frequency limit of an under¬ 
lying theory. Concluding that a system behaves hydro dynamically does not provide 
any information about the underlying physical processes. Also, for some cases, it is 
debatable if the model is initial- or final-state. With that being said, we will adopt 
the usual classification and identify models as non-hydro dynamic (initial-state) or 
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Fig. 18. Space-time diagram showing the causal relation between particles a and b. Particle a and 
b are causally connected to the red and blue patches respectively. The correlation between them 
must have formed in their overlap at a proper-time before Tq . 


hydrodynamic (final-state) in order to classify different approaches conveniently. 

All theoretical models that we consider have a common theme that makes them 
particularly interesting; correlations between pairs of rapidity-separated particles 
must have formed at early times.^^ In particular, if two correlated particles freeze 
out at a proper time ^ then causality dictates that the correlation must have 
formed at an earlier proper time constrained by 




(6) 


where i/a^b is the momentum-space rapidity of particles a, b. This relation can be 
derived by considering two particles at space-time rapidities rja and 775 . Assume 
particle a emits a signal at time ta = Tq cosh rja which is subsequently received 
by particle b at time 4 = q cosh 775 . The distance the signal must travel is 
ka — ^b\ = I^F.o. sinh 775 — Tq sinh77a|. Causality restricts \zb — Za\ < c{tb — ta) and 
therefore ^ exp - ^a|)- 

The factor of 1/2 in equation is obtained by introducing a third mutually 
correlated particle at rapidity 77 c. If this third particle were to emit two signals at 
time Tq particle a and b can be correlated as long as ^ exp {—\r]a — Vc\) ctnd 

Tq < Tp Q exp {—\r]b — Vc\)- We arrive at the most stringent condition for when 
the third particle, particle c, is located at mid-rapidity, {r]a + Vb)l‘^^ and find 




(7) 


The final step in arriving at the desired result of equation is equating the 
space-time rapidity with the momentum-space rapidity (z.e. 77 ^ 77 in the notation 
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Fig. 19. Yano-Koonin source velocity, Yyx versus pair momentum Ytttt demonstrating approximate 
boost-invariant longitudinal expansion.!^ 


used throughout this section). This is natural at high energies; Lorentz contraction 
makes the nuclei infinitely thin in the center of mass frame and there is no dimen¬ 
sionful scale in the longitudinal direction. Under a boost-invariant expansion the 
longitudinal flow of the system has the scale invariant form vl = zjt and therefore 
r] = y. 

However, it is not a priori the case that this scaling is maintained by the subse¬ 
quent dynamics. Once an additional scale is introduced the strict scaling relation 
Vl = zjt may no longer hold. For example, in thermal equilibrium the particle 
distribution has the form 

/ ~ exp (^-^ cosh(y - 77 )) (8) 

A thermal smearing of about one unit in rapidity can be expected as long as 
m±lT > 1. On the other hand, if pr T the correspondence between r] and 
y is completely lost. 

The best empirical evidence for a rigid correlation between the space-time ra¬ 
pidity 77 and momentum-space rapidity y comes from studying the Yano-Koonin- 
Podgoretskh (YKP) parameterizatioiSEH of measured Hanbury Brown-Twiss 
(HBT) correlation functions. The relationship between the Yano-Koonin source 
velocity, Yyk = tanh“^uvx, and pair momentum ^ log ^ 

allows for testing the longitudinal expansion.!^ For a non-expanding source Yyk 
would be independent of Y^^^. For a source undergoing a boost invariant expansion 
Yyk = and this appears to be satisfied by the available data at RHIC as 

shown in figure [T^ 

There seems to be a compelling case that the bulk of particle production is 
approximately boost invariant and the above causality arguments are well founded. 
The proceeding discussion has far-reaching implications on our understanding of 
long-range correlations. 
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The PHOBOS collaboration has observed that the near-side ridge signal in 
heavy-ion collisions persists up to six units in rapidity separationEHUll Using a 
very conservative estimate of ^ ^ 10 fm/c these correlations must have been 
established well before r ^ 0.5 fm/c. Long-range rapidity correlations thus probe 
early space-time dynamics and are therefore particularly well suited for studying 
novel QCD processes. For example, it has been argued that long range correlations 
may be a useful tool for probing high energy evolution of multi-parton correlations 
in the nuclear wavefunction.^^^ 

4.1. Final-state interactions: hydrodynamics in small systems 

From the earliest days of RHIC, hydrodynamics has been very successful in describ¬ 
ing the bulk properties of heavy-ion collisions Ideal hydrodynamics, which 
neglects viscous corrections, was able to reproduce the trends seen in experimental 
data, in particular the centrality and transverse momentum dependence of particle 
spectra and elliptic flow. Subsequent developments in the description of the initial 
state, in particular the inclusion of event-by-event fluctuations, and in the extension 
of the formalism to include off-equilibrium corrections like those introduced by a 
non-zero shear viscosity, have advanced the theoretical modeling significantlylHlll231 
Quantitative agreement of viscous relativistic fluid dynamics with a wide range of 
experimental data is excellent (see Fig.[2Q|), and recent improvements such as the 
inclusion of bulk viscosity and microscopic transport for the later stages of the 
collision improve the agreement even more.^^^ 



Fig. 20. Model calculations using a constant shear viscosity to entropy density ratio p/s = 0.18 
for Vn vs. colli sion cen trality in Pb-Pb collisions at UUUv = 2.76 TeV compared to experimental 
measurements 112£J112I1 


The similarity of experimental data for multi-particle correlations in small sys- 
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terns with those in heavy ion collisions (see Sec.[^ suggests that the physical pro¬ 
cesses responsible for the observed correlations may be the same in both systems. 
More explicitly, hydrodynamics could also be the appropriate model to describe 
high multiplicity p-p and p-A collisions. 

While the applicability of hydrodynamics in very small systems should definitely 
be questioned (and we will discuss this issue in detail below), the earliest discus¬ 
sion of (ideal) relativistic hydrodynamics to describe nuclear collisions involved p-p 
collisions Following the early argument by Landau, hydrodynamics should 
be applicable if at the moment of collision, a large number of “particles” is created, 
the “collision mean free path” in the created system is small compared with its 
dimensions, and a statistical equilibrium is established. We use quotation marks 
because the definition of particle number and mean free path may be ambiguous in 
a strongly interacting system. 

More recently hydrodynamic calculations aiming at the prediction and descrip¬ 
tion of experimental data for small collision systems at RHIC and LHC have become 
availableIn particular, p-A and d-Au collisions at high energy were 
explored early on, and predictions showed qualitative agreement with subsequently 
published experimental data (albeit at different collision energies). 

In the following we will discuss the ingredients of these calculations, emphasizing 
the differences between the calculations of various groups. We will see in particular 
that results for anisotropic flow are highly sensitive to assumptions about the initial 
state and its fluctuating structure. 

4.1.1. Modeling of the initial state 

The most important ingredient in the hydrodynamic description of small collision 
systems is the modeling of the initial geometry. Fluctuations in the initial energy 
density that deviate from sphericity propagate hydro dynamically and generate non¬ 
zero flow harmonics. In practice the initial state geometry is generated by doing 
event-by-event calculations where fluctuations of density distributions inside the 
incoming nuclei or hadrons are taken into account. In many calculations these fluc¬ 
tuations are dominated (or exclusively given) by the random positions of nucleons 
in a nucleus. However, in collisions of a single proton with a large nucleus the 
substructure of the proton may become important. 

We present the details of the initial state model calculations that have been 
employed for the description of p-A collisions below. These are the Monte Carlo 
(MC) Glauber model,II^the IP-Glasma model jMEIll and EPOS.I^^ 

Monte Carlo Glauber model We follow the discussion in the literature to in¬ 
troduce the implementation of the Monte Carlo (MC) Glauber model for p-A col¬ 
lisions 

In the MC Glauber model, which has been fairly successful in describing ex¬ 
perimental data from heavy ion collisionsthe proton and the energy deposited 
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per wounded (participating) nucleon are assumed to be spherical. Hence, in p-A 
collisions the fluctuating shape of the interaction region is solely due to the random 
positions of struck nucleons in the heavy nucleus. 

The specific implementation we discuss herd^^ uses GLISSANDOEHwhere the 
entropy deposited per wounded nucleon (or binary collision) fluctuates according 
to a gamma distribution in order that the final multiplicity distribution reproduces 
the empirically observed negative binomial. The width of each Gaussian shaped (in 
transverse space) contribution to the entropy density is chosen to be 0.4 fm. 

To determine whether a nucleon was wounded, a Gaussian wounding profile is 
used. This means that a nucleon undergoes a collision with probability 

P{h) = A exp(—7rA6^/crNN), (9) 


where b is the distance to the nucleon it is potentially colliding with, <jnn is the 
nucleon nucleon cross section, and A is a constant chosen to be 0.92. 

This method leads to a fairly large interaction region, much larger than the size 
of the incoming proton. This is illustrated in Fig. 21 a). If one chooses the alternative 
method of depositing energy around the interaction point of two nucleons, the 
average size is reduced by approximately 40 % 

In 3+1 dimensional simulations pHEZI the initial entropy density needs to be 
defined in the space-time rapidity direction as well. Typically one assumes the fac¬ 
torization of the transverse and longitudinal distributions. Assuming that the trans¬ 
verse distribution is independent of rapidity within some range around mid-rapidity 
is important to reproduce the experimentally observed ridge structure: Having the 
same initial transverse geometry at different space-time rapidities, leads to corre¬ 
lations between the flow pattern of particles emerging with different momentum 
rapidities. The typical profile in space-time rapidity is given by a combination of 
a linearly dropping central region and two half Gaussians at large (positive and 
negative) rapidities: 


p(7?s) = exp -%)) (1 ± 6>(2/beam ±??s) • (10) 

The linear drop away from the nucleon introduces a torque effect and some small 
fluctuation in the longitudinal direction. The parameters typically choserlSlEZI in 
calculations for p-A collisions at ^/s = 5.02 TeV are 770 = 2.5, = 1.4, and 

^beam = 8.58. Another parameter is the initial time that the hydrodynamic calcu¬ 
lation starts, typically chosen between tq = 0.2fm/c and tq = 0.6fm/c. 

Since in the hydrodynamic framework the final flow anisotropy is driven by the 
initial state eccentricity, we define it here and will refer to it later in the text when 
discussing how different initial state models generate different geometries. The n’th 
order eccentricity is given b^JSHIlll 


^Also using the black disk approximation, where a collision only occurs iib < R = leads 

to smaller interaction regions. 
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Models of initial states in pA collisions 


Glauber-like 



a) 


IP-glasma 



eccentric proton pomerons 



0 d) 


Fig. 21. Typical geometric configurations generated by different initial state models for p-A 
collisions: a) MC-Glauber with wounded nucleons as sources b) IP-Glasma with round proton c) 
IP-Glasma with an eccentric proton d) Pomerons. 


^n 


r” cos(n(/)))^ + (r" sin(n(/)))^ 

(F) 


( 11 ) 


where (•) is the energy density (or entropy density) weighted average, and r and (j) 
are the polar coordinates in the plane transverse to the collision axis. 

IP-Glasma model The Impact Parameter (IP) dependent Glasma modeEEHlis 
based on the color glass condensate frameworlPI^ and uses the IP-Sat modeP^to 
constrain the impact parameter and gluon momentum fraction dependence of the 
dipole cross section. This, together with fluctuating nucleon positions for heavy ions 
determines the initial distribution of color charges, which then enter the currents 
in the Yang-Mills equations that determine the initial gluon fields. The assumption 
in this model being that at high energies the QCD coupling constant is smal|^and 
the occupation numbers of gluons are non-perturbatively large (^ l/o^s)- In this 
limit the dynamics can be approximated classically by the Yang-Mills equations 
describing both the initial configuration of gluon fields and their dynamical evolu- 
tionJUmUl] heavy ion collisions at high energies these assumptions are likely 
to be valid. For smaller systems, such as p-A collisions, the approximations might 
be appropriate for high multiplicity events. Details on the implementation in A-A 
collisiondli^ll2II and p-A collision can be found in the literature 11^11211 Here we list 
the main ingredients of the IP-Glasma model: 


(1) IP-Sat parametrization of the saturation scale (5s(x,xt) and fixed parameters 

from Dis 


^The coupling is determined at the dynamically generated saturation scale Qs, which grows with 
increasing energy. 
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(2) Monte Carlo sampling of nucleon positions in the incoming nucleons and the 
impact parameter b. 

(3) Determination of the thickness function as a function of x^, T(xt), in each 
nucleus using information from (2). Using (1), this results in the distribution 
of Qs(xt) at the collision energy of interest. 

(4) Using a proportionality constant of order 1 that allows the conversion from Qs 
to the color charge density g‘^ii leads to the color charge density distribution 
in each nucleus. This proportionality constant can fluctuate to account for the 
effect of gluon number fluctuations. See discussion below. 

(5) Monte Carlo sampling of color charges (with color index a) for each nucleus 
from the color charge density distribution, assuming 

{p°‘{yiT)p^{yT)) = - yr) • (12) 

(6) Solution of the Yang Mills equation 

= J", (13) 

for each nucleus with the color currents given by 

= 5'^±p“(a;=F,XT)r, (14) 

where + and — are for right moving and left moving nuclei, respectively, are 
the generators of SU(3). 


(7) Numerical solution for the gluon fields after the collision at time r = 0+E2I1E121 

Ai — Ai 1 Ai 

(15) 

II 

^\a)A\b) 

, (16) 

O 

II 


(17) 

drA^ = 0, 


(18) 


where (A) and (B) denote fields of the two colliding nuclei. 

(8) Solution of the source-free Yang-Mills equation to obtain the time evolution of 
the produced gluon fields. 

(9) Landau matching of the gluon-field energy momentum tensor to obtain the 
relevant quantities for fluid dynamic calculations: energy density ^(xt) and 
flow velocities 14 ^(xt). 

Note that initial conditions obtained with this model are boost invariant, which is a 
good approximation for high energy collisions and rapidities not too far away from 
mid-rapidity. 

We will now discuss specific features of the model relevant for small collision 
systems. Due to the fluctuating color charges, gluon multiplicity distributions in 
the IP-Glasma model already include negative binomial fluctuationsE4lll2ZI It was 
however discovered that for p-p and p-A collisions additional fluctuations (as com¬ 
pared to the original modeP^ of the color charge density were needed to describe 
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the experimental multiplicity distributions (see also Fig. 35). These fluctuations can 
be understood as fluctuations of the gluon number for a given value of the satura¬ 
tion scale Qs. The comparison of gluon with hadron distributions is of course only 
qualitative and additional fluctuations from the hadronization mechanism are also 
expected to modify the result. Nevertheless, it is important to include the gluon 
number fluctuations to allow for rare fluctuations that generate the high multiplicity 
events of interest. 

Turning to the initial geometry, the IP-Glasma model generates interaction re¬ 
gions in the transverse plane that are dominated by the smaller of the two colliding 
objects. This is because energy is only deposited where both nuclei’s gluon fields 
are non-zero (otherwise the fields are pure gauge fields and indistinguishable from 
the vacuum). For a proton colliding with a large nucleus (not too close to its edge), 
the interaction region is thus strongly affected by the shape of the proton. This be¬ 
havior leads to much smaller interaction regions than in the Monte Carlo Glauber 
model as illustrated in Fig. 21 b). Furthermore, if the proton is assumed to be round 


(up to fluctuations of the color charges p^), initial eccentricities as defined in (11) 
are very small. 


EPOS The EPOS modeE3E211E23 has successfully been applied to describe exper¬ 
imental data from both p-p and p-Pb collisions. EPOS includes both a fluctuating 
initial state and hydrodynamic evolution, as well as jets and the interaction between 
jets and the medium. The initial state calculation is based on multiple scatterings, 
where initial state radiation is included. Thus, the model deals with parton lad¬ 
ders or “Pomerons”. The Pomerons, whose number determines the multiplicity, are 
placed in between two colliding partons, which are distributed around the center of 
nucleons according to an exponential distribution. In addition, the pomerons carry 
a transverse momentum (generated by the initial scattering) and thus move trans¬ 
versely even before the coupling to hydrodynamics occurs. EPOS does not assume 
boost invariance but produces a smoothly changing longitudinal energy density 
distribution. It should be noted that EPOS also needs to include saturation effects 
in order to not violate unitarity. It is done by introducing an energy dependent 
saturation scale for each Pomeron.l^ 


Holographic pomerons Eollowing the AdS/CET conjectur^^^ a reggeized and 
unitarized scattering amplitude was derived for A/' = 4 SYM in the limit of large 
N, strong coupling Qym^ ^ ^ high energy Iii^ it was recognized that a 
single closed string exchange contribution to the eikonalized dipole-dipole scattering 
amplitude yields a Regge behavior of the elastic amplitude, with features consistent 
with experimental data. 

This Regge behavior is driven by worldsheet instantons, which describe the 
Schwinger mechanism for string pair creation by an electric field, where the longi¬ 
tudinal electric field at the origin of this non-perturbative mechanism is induced by 
the relative rapidity of the scattering dipoles. 
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This stringy Schwinger mechanism can generate quantum fluctuations of virtual 
strings, which take the form of thermal fluctuations with an effective temperature 
related to the inverse string circumference and therefore impact parameter; at the 
center of the string T~^ ~ 27r|bx|. If the impact parameter is small enough there 
may be a region in the middle of the string where the effective temperature is 
above the deconfinement temperature suggesting that the string description should 
be replaced by a plasma phaseJ^^ 

It has been argued that near this critical impact parameter the increased quan¬ 
tum fluctuations reduce the string tension and self-interactions among the strings 
cause an implosion and formation of a string ball^^ A more quantitative study of 
string self-interactions and string ball formation was explored inlHllHD At effective 
temperatures larger than Tc the string ball may form a strongly coupled high density 
plasma that expands hydrodynamically. One of the key features of the holographic 
pomerons that may distinguish them from other pictures of plasma formation is the 
sudden onset of hydrodynamic behavior at a critical impact parameter. However, 
no explicit hydrodynamic calculations have been performed using this framework. 


4.1.2. Viscous hydrodynamics and its limits of application 

Hydrodynamics is an effective theory for the long wavelength modes of a system 
able to describe the interactions of the low momentum part of the particle spec¬ 
trum if, from a kinetic theory point of view, the mean free path of the particles 
is significantly smaller than the system size. This lead many to criticize the use 
of viscous fluid dynamics in small colliding systems and indeed one should always 
check whether hydrodynamics is being applied within its domain of validity. In this 
section we will first introduce the basic ingredients of viscous relativistic hydrody¬ 
namics and quantify the domain of validity more explicitly. 

The viscous relativistic fluid dynamic equations for small systems are identical 
to those employed to describe heavy ion collisions. Various recent reviews cover the 
details of both the equations and their derivation, as well as the implementation in 

Here we summarize the most important aspects. Relativistic viscous hydrody¬ 
namics is formulated as a gradient expansion to take into account deviations from 
local thermal equilibrium or ideal hydrodynamics. Due to the acausal behavior of 
the relativistic equations expanded to first order in gradients, numerical implemen¬ 
tations require the inclusion of terms at least of second order in the gradients. 
Commonly used are Israel-Stewart type equat ionilHllli21 that were derived using 
the 14 moment method and the second moment of the Boltzmann equation to 
determine the hydrodynamic equations for the dissipative currents. Improvements 
that provide a more systematic expansion were derived more recently.f^^ 

Many simulations concentrate on shear viscosity, but some also include bulk vis¬ 
cosity, and possibly the coupling between shear and bulk terms. In fact, the effect of 
bulk viscosity has been found to be larg^ISH for certain observables, including spec- 
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Fig. 22 . Space-time evolution of the Knudsen numbers in p-Pb collision at the LHC, with r]/s = 
0.08 in the QGP phase and a te mperature dependent 77 /s in the hadron gas phase, and dN^h/dr] = 
270. (a) Kne and (b) Kne^ 


tra and anisotropic flow. However, a large uncertainty arises from off-equilibrium 
corrections to the particle spectra when converting hydrodynamic quantities to 
particles at freeze-out, in particular when including bulk viscous corrections. 

Large corrections to the equilibrium distribution functions (typically appearing 
at transverse momenta Pt ^ 2 GeV) are one indication that the fluid dynamic 
approximation we are using breaks down. Another, more systematic quantification 
of the validity of viscous fluid dynamics is provided by the study of Knudsen and 
Reynolds numbers. The Knudsen number is defined as the ratio of a microscopic 
to a macroscopic scale in the system. Various choices for these scales are possible, 
and a variety of them has been studied in both heavy ion collisions and smaller 
systems.fi^ For example, the macroscopic scale can be defined as the inverse of 
the expansion rate with the flow velocity, while the 

microscopic scale can be described by the shear relaxation time, /micro = which 
in a dilute gas is proportional to the mean free path, ^ Amfp- So, Kne = 
or alternatively we can define with where eo is the 

energy density. 

For hydrodynamics to be valid, the Knudsen number should stay significantly 
below one at all points in space time. Whether this condition is fulfilled was tested 
for various choices for Lmacro^^^and we show examples for p-A collisions in Fig. 

One can see that around T = 100 MeV the fluid dynamical description is no 
longer applicable, with all Kuq values being above 1.5. The results shown are for a 
constant rj/s = 0.08 in the QGP phase - if a temperature dependent r]/s was used 
in the QGP phase, the situation would worsen and fluid dynamics would be out 
of its domain of applicability even in the early stages of the evolution. While large 
Knudsen numbers indicate the breakdown of hydrodynamics it does not preclude 
large final state effects. However, it does mean that another framework outside of 
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hydrodynamics might be more appropriate to model the system when the Knudsen 
number is large. 

4.1.3. Results from hydrodynamics in small systems 

We now discuss the results obtained for the Fourier harmonics in small collision 
systems within various model calculations, beginning with the existing predictions 
for proton-proton collisions at the LHC and continuing with proton/deuteron on 
heavy ion collisions. We will then go into more detail with the discussion of more 
recent, very sophisticated calculations and their results for mean transverse mo¬ 
mentum, Fourier harmonics of charged hadrons and identified particles, as well as 
HBT radii. 

Proton-proton collisions Wide ranges of predictions for elliptic flow in proton- 
proton collisions were published early on - in particular in expectation of LHC 
results at 7 and 14 GeV. We list several calculations and predicted values for the 
elliptic flow in p-p collisions in the following. 

Ideal hydrodynamic calculationindicated that V 2 ^ 0.035 in p-p collisions, 
and that any value above 0.02 in minimum bias collisions would indicate an ex¬ 
tremely small viscosity (below the KSS bound of l/47r) of the medium. 

In a string percolation model that includes a directional dependence of the 
sources, a sizable V 2 was predicted for p-p collisions .^23 Various scenarios including 
different shapes of the proton have been studied^ and estimates of V 2 for p-p 
collisions depended strongly on the assumed shape, ranging from zero to almost 
1.5% for the more realistic models. Similarly strong dependencies on the proton 
profile were found in hydrodynamic simulations,^^^ where very large values for the 
minimum bias V 2 ^ ranging from 6% to 25% are quoted. 

High multiplicity events were predicted to have an integrated V 2 of around 10%, 
but again with a wide spread depending on the initial state configuration, when 
using hydrodynamics after initializing with flux tubes obtained in a constituent 
quark model.li^ A similar model using partonic interactions determines eccentric¬ 
ities compatible with V 2 values ranging from 10% to 20%, depending on model 
parameters, when assuming a given relation between eccentricities and elliptic flow 
values determined from hydrodynamic calculations.^^^ 

Mimicking quantum fluctuations by introducing an eccentricity to the inter¬ 
action region, and then performing hydrodynamic simulations, a V 2 value of 3% 
was found if the induced eccentricity has average value of about 17% in another 
study. 

Using DIPSY, a Monte Carlo event generator based on the QCD dipole model, 
to generate the initial conditions and a fixed relation between V 2 and eccentricity, a 
value of 6% for V 2 in p-p collisions at LHC was estimated, being fairly independent 
of the multiplicity.^^^ 

In the EPOS framework that performs hydrodynamic expansion based on flux- 
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Fig. 23. Predictions for the elliptic flow coefficient V 2 from a 3+lD viscous relativistic hydrody¬ 
namic simulation for p-Pb (at ^/s = 4. 4 TeV) a nd d+Pb (at ^/s = 3.11 TeV) collisions, compared 
to ALICE results for Pb+Pb collisions 

tube initial conditions as described above, both the experimentally determined 
mean transverse momentum and measures of the source size were reproduced?^ 
Furthermore, the ridge structure was reproduced within this model?!^ This in¬ 
dicates that the elliptic flow component is correctly described by EPOS in p-p 
collisions. 

The conclusion to be drawn from this extremely wide range of predictions is 
the strong model dependence, in particular on the initial state for hydrodynamic 
calculations, that one faces when computing V 2 in proton-proton collisions. 


Proton-heavy ion collisions Early predictions for flow in d-Pb and p-Pb col¬ 
lisions from hydrodynamic simulations were also present ed.^^^ In this work MC- 
Glauber initial conditions, as described in Section [4.1.1[ were used together with 
3+1 dimensional viscous hydrodynamic simulations to obtain results for particle 
spectra and flow harmonics as functions of transverse momentum and pseudo ra¬ 
pidity. 

The elliptic flow is found to be 3%-4% in p-Pb collisions, with little centrality 
dependence (Eig.|^, significantly smaller than for Pb-Pb collisions with the same 
multiplicity. Eor the d-Pb system, the elliptic flow is significantly larger, increasing 
for central collisions, and reaching almost 10%. 

The ridge structure in two-particle correlations was observed in similar hydrody¬ 
namic calculations for p-Pb collisions at ^/s = 5.02 TeV.Eor the largest rapidity 
gaps considered, V 2 integrated over 0.1 GeV < Pt < 2 GeV was found to be 4%, 
while the corresponding was about half as big. 
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As discussed before, in these calculations the initial interaction region turns out 
to be rather large when using scaling of the initial energy or entropy density with 
the participating nucleons. The size and lifetime of the collective source formed in 
central p-Pb collisions is found to be 3 — 4 fm. 

More quantitative results on the size of the emission source using HBT tech¬ 
niques were subsequently present ed.l^^ Here, an alternative model for the initial 
distribution was also considered. This “compact” variant locates the sources in the 
center of mass of nucleon-nucleon collisions. The average rms radii for the two sets 
of initial conditions considered were 1.5 and 0.9 fm, respectively. 

Results obtained for HBT radii as a function of multiplicity were presented 
in comparison to heavy ion and p-p data.^^^ It was found that the p-Pb results 
for Rside, Rout^ aud Riong were much closer to the heavy ion results at a given 
multiplicity than to the p-p results. Rout is reduced by 25% when using the compact 
initialization, while the other radii are little affected. We note that the difference in 
the initial size of the system has a large effect on the final measured HBT radii. In 
particular, it was showrP^ that the IP-Glasma model (see Section 4.1.1) produces 
significantly smaller initial energy density distributions in the transverse plane than 
either of the MC Glauber models discussed above. In particular, initial sizes in p- 
Pb collisions were demonstrated to be much closer to those in p-p collisions than 
Pb-Pb collisions.One would expect that this leads to similar differences for HBT 
radii obtained in hydro calculations using the different initial state models. 

Very detailed analyses of flow in p-Pb collisions using the MC Glauber model 
including negative binomial fluctuations, which are necessary to reproduce the mul¬ 
tiplicity distribution, and in particular its tail at high multiplicity, have been carried 
out following the first exploratory works.^^^ Here, the authors employ hydrody¬ 
namic simulations followed by statistical freeze-out. The negative binomial fluctua¬ 
tions in the initial state, which cause the deposited entropy per participant nucleon 
to fluctuate, turn out to be necessary in order to describe the experimental data 
by increasing the initial eccentricities and thus the resulting U2{2} and U2{4}. The 
successful description of ATLAS datgP^ in this case is presented in Fig.[^ 

Gomparing to predictions discussed above and shown in Fig.|^ we conclude 
that (differences from the employed method to determine V 2 aside), large differences 
result from using different initial state models. 

More detailed analyses have demonstrated that various features of the data are 
naturally reproduced by hydrodynamic calculations, even on a quantitative level. 
For example, the mass splitting of the mean transverse moment unCUlMI and the 
V 2 of identified particleJnUZH in p-Pb collisions was presented using different initial 
state descriptions. 

We note that hadronic cascade models, such as UrQMD, also reproduce the 
characteristic mass splitting. UrQMD could however not describe the centrality 
dependence of the elliptic flow or get a real valued U2{4}.II^The latter is naturally 
achieved in hydrodynamic calculations, where all higher order cumulants (from 4, 
6, 8 particles, etc.) are equal since all particles are correlated with one another. 
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Fig. 24. The elliptic flow coefficient of charged particles for \r]\ < 2.5, 0.3 GeV < pr < 5.0 GeV 
from the cumulant method and i;2{4}, and from the di-hadron cor relat ion function measured 

by the ATLAS Gollaboration,compared to hydrodynamic calculations.^^ Glauber+NB includes 
negative binomial fluctuations in the initial entropy density distribution. 


The IP-Glasma model produces an initial transverse energy distribution having 
the shape of the assumed proton shape, and early calculations using a spherical 
proton therefore yielded a very small V 2 in p-A collisionsHowever, the pro¬ 
cess of energy deposition is less ambiguous as in the MC Glauber model, where 
very different results can be achieved by using different prescriptions for the initial 
energy/entropy density distributions. 

Within the IP-Glasma + hydrodynamics description it is thus clear that the 
results for azimuthal anisotropies in p-Pb collisions are sensitive to the shape of the 
proton and its fluctuations event by event 

In conclusion, all aspects of the experimental data are at least qualitatively 
and for the most part also quantitatively described by hydrodynamic calculations. 
Unfortunately, in addition to the unknown transport parameters, results are highly 
sensitive to the assumptions made in modeling the initial state. Furthermore, one 
must scrutinize the applicability of hydrodynamics in these very small systems due 
to large Knudsen numbers throughout the evolution. 

Some support for the applicability of hydrodynamics comes from calculations 
at large coupling, where it was shown that viscous hydrodynamics describes early 
on the evolution obtained using AdS/GFT correspondence in a system of colliding 
shock waves. 
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System sizes in pp, pA and AA collisions 
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Fig. 2 5. T he experimental radius i^inv in P-P? P-Pb? cind Pb-Pb collisions from the ALI CE collab- 
oratioi®^ (left) compared to the radius rmax in the IP-Glasma (+Music) model.EMl The lower 
end of the band indicates the initial size, the upper end the maximal size reached during hydro- 
dynamic expansion. Initial sizes are similar in p-p and p-A collisions, while A+A collision systems 
at the same multiplicity are generally larger. This is a trend reflected also in the experimental 
HBT radii shown on the left. 


HBT radii from hydrodynamics As for the results on the azimuthal anisotropy, 
results for the system size, mainly characterized by HBT radii (discussed in Section 


3.5), vary largely depending on the initial state model. While the IP-Glasma model 


produces initial sizes that are similar in p-p and p-A collisions, but large in A-A for 
the same multiplicity, the MC-Glauber model produces larger p-A systems, closer to 
those in A-A collisionsas illustrated in Fig. 21 While HBT calculations within 
the IP-Glasma model have yet to be performed, the general trends of the initial size 
of the system and the measured HBT radii agree very well. In Fig.we compare 
the experimental data on the left to results for rmax on the right. We define rmax as 
the (angle averaged) radius where the system reaches a minimal threshold energy 
density Smin = This radius by definition depends on the choice of Smm- This 

choice however only affects the overall normalization of rmax; h does not affect the 
dependence of rmax on the number of charged particles A^ch- There is also some 
dependence on the infrared cutoff m, which has been chosen to be on the order 
of Aqcd- This dependence can mostly be absorbed into an overall normalization 
constant.^^^ 

In summary, the system sizes obtained from experimental HBT measurements 
are generally compatible with the fluid dynamic picture, while the details, once 
again, depend strongly on the initial state and its size. 


Deuteron/^He-heavy ion collisions Small systems with drastically different 
geometries from p-A collisions could help distinguish whether the origin of the 
azimuthal anisotropies measured in p-A collisions could indeed stem from the hy¬ 
drodynamic response to the initial geometry. 

Therefore, the PHENIX experiment has measured V 2 in d-Au and V 2 and 
in ^He-Au collisions. Hydrodynamic calculation£lllElHl predicted a somewhat en- 
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Fig. 26. Experimental results for V 2 and V 3 in ^He+Au from PHENI}<P^ compare d to three differ¬ 
ent hydrodynamic calculation s: A prediction using MC Glauber initi al con ditions'^ results from 
superSonic including pre-flow,and IP-Glasma+MusiC calculation,^^ with increased viscosity 
of 77 /s = 0.18 compared to the original prediction. 


hanced V 3 for the ^He-Au system because of the enhancement of significantly tri¬ 
angular events. Indeed, hydrodynamics can describe the v^s also in these two sys¬ 
tems EM1EIZ1E13 While some predict ion had to be amended with a “pre-flow” 

component, other^^^ need an adjustment of the shear viscosity. Both calculations 
then lead to good agreement with the experimental dataEUEHl In Fig. 26 we show 
a compilation of results for V 2 and vs in ^He-Au collisions. While some of the cal¬ 
culations, which are predictions, can be improved by varying the shear viscosity, 
general agreement with the experimental data is very good, for the very different 
initial conditions. superSonic and IP-Glasma+MusiC seem to have a somewhat too 
large ratio of ^ 2 /^ 3 . It would be desirable to also experimentally measure vs in d- 
A and p-A collisions to see whether it is indeed enhanced in ^He-A collisions as 
expected from the initial geometry and confirmed in hydrodynamic calculations.^^^ 


4.1.4. Parton transport models 

Parton transport models describe final state interactions microscopically, including 
various kinds of processes (elastic, inelastic) and different implementations of cross 
sections (isotropic, pQCD, ...). If these cross sections are large enough, the micro¬ 
scopic description should approach the hydrodynamic limit, and we expect similar 
results for flow observables as in (viscous) hydrodynamic simulations. 
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Fig. 27. The transverse momentum dependence of the elliptic, ^ 2 , and triangular, v^^ anisotropy 
coefficients in p-Pb (upper panel) and Pb-Pb collisions (lower panel) from the AMPT mode l (op en 
symbols) compared to experimental data from CMS (filled symbols). Figure adapted from.E^ 


AMPT The AMPT modePSlis based on HIJINGP^ to generate the initial state 
and the parton cascade ZPCP^ that simulates the final state interactions. One 
important ingredient that allows AMPT to produce enough final state collectivity 
is the so called string melting mechanism, which converts all initial minijets and 
soft strings into quarks and anti-quarks. These undergo elastic scatterings with a 
partonic cross-section which is controlled by the strong coupling constant and the 
Debye screening mass. After hadrons are formed using a coalescence model, they 
also interact further. 

p-Pb (and p-p) collisions have been studied within the AMPT model inll551llMl 
It was found that with a modest elastic cross section of cr = 1.5 — 3 mb long range 
two-particle correlations were well described and the computed V 2 and agreed 
well with the experimental data from p-Pb collisions (while V 2 in Pb-Pb collisions 
was underestimated by the model) - see Fig.[^ 

A particularly interesting observation of this stud^ff^ was that on average each 
parton undergoes only two collisions (and the hadronic cascade had a negligible 
effect on the results). It is surprising that such a small amount of interactions could 
lead to the build-up of collective flow. 

A possible explanation of the physical process in the AMPT model was given in 
a recent work,^^ where the details of what partons carry the anisotropy and how 
it depends on the number of their scatterings were analyzed. It was found that an 
escape mechanism seems to be responsible for the produced anisotropy in AMPT. 
The main ingredient of this is the initial spatial anisotropy of the medium and the 
resulting greater probability for partons to escape along the direction in which the 
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medium is shorter. 

This effect is drastically different from the build-up of hydrodynamic flow, and 
if it is the main effect, it is somewhat surprising that the transverse momentum 
dependence of the Vn coefficients is so well reproduced by AMPT. 

4.2. Conclusions on final state interactions in small systems 

In conclusion we can state that it is possible to describe all characteristic features 
measured in p-p and p/d/^He-A collisions with models based on the collective 
response to an initial state geometry. In particular hydrodynamic models can re¬ 
produce the azimuthal anisotropies of charged hadrons Vn^ the mass splitting of 
the mean transverse momentum and V 2 for identified particles, and the HBT radii. 
However, for all observables there is a significant uncertainty related to the little 
knowledge we have about the initial state in small systems. Depending on the as¬ 
sumptions made about how the initial shape of the system is generated, final results 
can vary dramatically. Apart from this caveat, the Knudsen number can take on 
rather large values in viscous hydrodynamic simulations of small systems. Thus, the 
quantitative results could be plagued with errors due to running hydrodynamics at 
least partly outside its domain of applicability. 

4.3. Multi-particle production in the Color Glass Condensate 

Introduction to the Color Glass Condensate A first attempt at understand¬ 
ing the nature of multi-particle correlations should come from studying the pro¬ 
duction processes dictated by QCD. Traditionally this is done in the framework of 
collinear factorization where the strong interaction processes are factorized into a 
product of cross sections convoluted with parton distribution functions. This ap¬ 
proach is valid in the Bjorken (or short distance) limit, 

Aqcd Q ^ ^ (19) 

where Q is the kinematic energy scale of the scattering process, and x the Bjorken 
X, z.e., the longitudinal momentum fraction of the parton in the infinite momentum 
frame. Requiring the independence of an energy scale y separating of the long dis¬ 
tance physics encoded in the parton distribution functions from the short distance 
QCD matrix elements leads to the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi^JpUmini renormalization group equation. By construction the DGLAP equa¬ 
tion resums collinear logarithms of the form as In (Q^/Aqqj^). 

However, most of the particle production and the bulk of QGD phenomena 
occurs at low Q and therefore we consider QCD in the Regge-Gribov (or high 
energy) limit, 

Aqcd <^k± x^/s < (20) 

where another class of large logarithms, Og ln(l/x), must be resummed, in this case 
by the Balitski-Fadin-Kuraev-Lipatov (BFKL) equation UUllH^ The rapid rise in 
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Fig. 28. Left: Longitudinal momentum frac tion, x, times the unpolarized parton distribut ions 
obtained in NNLO NNPDF2.3 global analysi^^at scale /jp = 10 GeV^. Figure adopted from.^^ 
Right: Evolution of the saturation scale Qs and the schematic structure of a nucleus with along 
the horizontal axis and Bjorken x along the vertical axis. The increase of the gluon number at a 
fixed transverse size scale along the vertical axis eventually leads to the phenomenon of saturation. 


the gluon distribution xG{x^ Q‘^) at small x for a fixed (see figure , growing 
approximately as a power of x~^'^ can be understood within the BFKL framework. 
The rapid growth is moderated by the fact that gluons having transverse momenta 
less than a semi-hard saturation scale Qs{x) have a maximum occupation oil/ag. 
A parametric estimate for the onset of saturation can be made by considering 
the maximum number of partons that can be packed in a proton of size 
The condition that the transverse area occupied by a parton Sx/xG) is on 
the order of the gluon fusion cross-section a ^ olsIQ"^ defines a saturation scale 
Ql ^ as xG{x^ QD/aS'x- The argument can be extended to nuclei; A nucleus with 
atomic number A has transverse size S± ^ and its parton distribution scales 
approximately as xGnucieus ^ ^ * xG resulting in 

ry A A^/'^ 

Q%{x) ^-^xG{xX%) - ^ (21) 

thereby enhancing saturation effects by a factor as large as 6 in nuclear 

targets such as lead (A = 208). The evolution of the saturation scale as a function 
of and x is schematically shown in the phase diagram of high energy QCD in 
figure 

The Color Glass Condensate is an effective field theory description of the sat¬ 
urated gluons in the Regge-Gribov limit. More details can be found in the recent 
review articlc.^^^ The degrees of freedom consist of strong color sources ^ 1/g 
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of the nuclei and the resulting classical gauge fields. The strong color sources vary 
event-by-event and are governed by a gauge invariant probability distribution, W 
with expectation values of operators schematically written as 

{O) ^ I [Dp,] [Dp,] W[p,]W[p,]0 [p„p,] . (22) 

The requirement that physical observables be independent of the scale separating 
the sources from fields results in the JIMWLK renormalization group equation. 
There is no known analytic solution to the JIMWLK evolution equation and while 
numerical solutions are available their use in phenomenology remains lim- 
ited.l^^ There exists a kinematic window in x where the nuclear wave-function can 
be derived in closed form. In the McLerran-Venugopalan modeP^^®^ x must be 
large enough, aslog{l/x) 1, such that quantum corrections entering at leading 
logarithmic order can be neglected and x must be small enough, x <C such 

that the small-x partons couple coherently to the large-x sources. In this limit, the 
weight functional has the following Gaussian form 


W[x^ p] = exp 



p^(x^)p^(x^) 


(23) 


where p\ = g^A/{2S±) is the average color charge squared per unit area per unit 
color degree of freedom. While the MV model is a useful phenomenological tool at 
moderate values of x ^ 10“^, at smaller x, as probed by the kinematics at the LHC, 
quantum corrections at leading logarithmic order inl/x must be incorporated. The 
difficulties with the JIMWLK equation arise from its hierarchical structure {e.g. the 
evolution of the two-point function depends on the four-point function). This can 
be circumvented by working in the large Nc approximation (in addition to the 
already imposed large nucleus approximation). In this limit the expectation value 
of a four-point function can be recast as the product of correlators of two point 
functions and, at leading order in the JIMWLK equation is simplified to the 
LO Balitsky-Kovchego\ElllESll (BK) equation 


^ /C(r, n) [T(ri, F) + T(r 2 , V) - T(r, Y) - T(ri, F) T(r 2 , F)] ,(24) 


In the above expression r 2 = r — ri and T(r, V) is the quark-anti-quark dipole scat¬ 
tering amplitude from a nucleus. In the context of the BK equation, the kernel JC 
amounts to resumming all terms of the form (a^ ln{xo/x))^ arising at any order in 
perturbation theory. It is well known that running coupling corrections qualitatively 
modify the small x evolution and analytic results determining the scale of the cou¬ 
pling that enters the evolution kernel are knownEMEHl The form of these running 
coupling corrections amount to two pieces; a piece resembling a running coupling 
correction to the LO BK Kernel,!!^ and a “subtraction term” that introduces new 
structures. In what follows we use the “Balitsky prescription” (equation (25)) for 
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-Y 

x=Xoe 




Fig. 29. Left: Un-integrated gluon distribution in the adjoint representation, plotted as 
q;s0/(S'_l A^c), with initial condition {Y = 0) corresponding to a min. bias proton along with 
subsequent evolution by rcBK to {Y = 4,8,12). Right: The saturation scale defined as the 
location of the maximum of k^(f)iY^k±_) as a function of T = Inixo/x) for initial conditions 
Ql = 0.168,1.008, 3.696 GeV^. 


the evolution kernel, 


/CBai.(r ri,r2) 


as(r)Vc 

TT 


2^2 


rtr: 


1^2 


1 (c^sYi) T 

d Uahi) ) 


1 (c^siA) 

d U4d) 



, (25) 


as it was shown that the “subtraction term” is numerically less importantP^ and 
will be ignored in this work. Numerical result^^^ for the full NLO BK equatiorP^ 
have recently become available demonstrating a sensitive dependence on initial con¬ 
dition with results not necessarily positive definite calling for a better understanding 
of the NLO BK equation before it can be used for phenomenology. A recent analysis 
that resums double logarithms seems to improve the situation at NLO.^^^ 

Solution of the running coupling BK equations (rcBK) requires a suitable initial 
condition at large-x. This review will focus on the results using an MV model-like 
initial condition with anomalous dimension 7 , 


T(r, V = 0) = 1 — exp 



7 

In 


( —— 
\rAQCD 



(26) 


The value of 7 and initial saturation scale (in the fundamental representation), 
have been constrained by global fits to available small-x {x < 0 . 01 ) lepton- 
proton scattering data.^^^ This work will use the parameter set h’ fronP^ having 
parameters Aqcd = 0.241 MeV, 7 = 1.119, Q^q = 0.168 GeV^. 

At large Nc the dipole scattering ampliutde can be identified with the uninte¬ 
grated gluon distribution in the adjoint represent at iorP^ 

i>(y,k_L) = f dh [l-T{r,Y)f . (27) 


For completeness we also point out that at large the Wilson lines can be ex¬ 
panded in powers of the sources yielding a relation between the average squared 
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color charge density and unintegrated gluon distribution 


^Y,k^) = {N^ - 1) 




(28) 


The numerical solution of the rcBK equation is shown in the left plot of figure 
The solid black curve is the initial condition defined in equation (26) and its 


evolution to T = 4, 8,12 is shown as the dashed, dotted, dash-dotted curves re¬ 
spectively. As we will argue at length below, the solitonic structure is crucial to our 
understanding of the ridge phenomenon. The non-monotonic behavior of the gluon 
distribution kinematically constrains gluons to be produced in similar directions. 
The transverse momentum dependence of the ridge will be controlled to a large 
extent by the maximum gluon occupation at the relevant values of x probed in the 
processes. 

The right plot of figure shows the location of the maximum of the quantity 
k±^{k±) as a function of Y or equivalently x = xoe~^. The quantity k±^{k±) 
represents the number of gluons at a given k± and its maximum value can be 
identified with a saturation momentum though it should not be confused with the 


saturation scale that enters into the initial condition in equation (26). The latter is 
the initial saturation scale at x = xq = in the fundamental representation. 

The solid purple curve in the right plot of figure corresponds to the wavefunc- 
tion shown in the plot on the left {i.e. having initial condition = 0.168 GeV^ 
which we will refer to as “minimum bias”). The rcBK equation has no impact pa¬ 
rameter dependence and the fits to lepton-proton scattering data (from which the 
value of QIq = 0.168 GeV^ is obtained) assumes that the impact parameter de¬ 
pendence factorizes from the dipole amplitude. Since the median impact parameter 
dominates inclusive scattering we will use this initial condition for proton configu¬ 
rations prototypical of minimum bias. 

In order to mimic rare fock space configurations generated in high-multiplicity 
events the initial saturation scale will be adjusted to account for the increased 
multiplicity. The right plot of figure also shows the evolution of the saturation 
momentum for configurations involved in high multiplicity events. The initial condi¬ 
tion having = 1.008 GeV^ corresponding to 6x the value of a min. bias proton 
will be argued to be the relevant scale need to describe high-multiplicity events in 
p-p collisions at the LHC. The initial condition with = 3.696 GeV^ (22x the 
min. bias proton) is used to describe a hot-spot configuration in the nucleus probed 
in central p-Pb events at the LHC. This large value is expected as it takes into 
account the additional enhancement from nuclear coherence. 


4.3.1. Two-gluon production for px ^ Qs- Glasma Graphs 

In the presence of strong classical sources the naive power counting of Feynman 
diagrams can change drastically. We consider two mechanisms for the production 
of a pair of mid-rapidity hadrons. The first proceeds via the production of two 
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Table 1. Power counting for multiparticle production in QCD in the dilute (left column) and 
dense (right column) limit of color sources. The top row shows the diagram responsible for mini¬ 
jet production (hadrons produced primarily back-to-back in azimuthal angle). The bottom row 
shows the “glasma” diagram responsible for the near-side collimation that is long-range in rapidity. 
Each three-gluon vertex or coupling to a bare parton is proportional to g. Each gluon coupling to 
a blob contributes 1 /g. 


gluons from a single t-channel gluon exchange as shown in the upper left amplitude 
in table We refer to this diagram as the Jet Graph, since the gluon pairs that 
are produced in this channel are predominately back-to-back in azimuthal angle. 
We should stress that we are interested in mini-jets (ie. soft dihadron correlations) 
and not jets in the usual sense studied in collinear factorization. 

The second class of diagrams shown in the lower left of Fig.has two gluons 
produced from two t-channel gluon exchanges. As drawn here the diagram appears 
disconnected and one might conclude does not contribute to an intrinsic two-particle 
correlation. This is not the case because the averaging over the color sources intro¬ 
duces non-trivial connections within the seemingly disconnected diagram. 

In the dilute limit each vertex (both the three-gluon and the coupling to the 
valence parton) come with a single power of the strong coupling constant g. After 
squaring the jet amplitude the cross section scales as a^. The glasma diagram on 
the other hand has six such vertices and the cross-section is proportional to and 
is therefore suppressed by in the dilute limit of both projectile and target. 

In the high energy limit the power counting changes dramatically. By high 
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energy limit we mean that k± ^ ^/s such that there is sufficient small-x evolution 
that the gluon densities are non-perturbatively large. The high-energy limit is also 
reached for large nuclei or rare collisions of hot-spot configurations of the projectile 
and/or target. The two tagged gluons now couple to the high occupancy ^ ^/ols of 
larger-x sources. The effective coupling of a gluon to the valence region is now 1/g 
instead of g. This is depicted in the right diagrams of table where each connection 
to a blob couples with 1/g. In the dense-dense limit the Jet Graph goes as while 
the glasma graph scales as aj‘^. For high multiplicity collisions one therefore expects 
the glasma graph as the main contributor to two-particle production. 

These two classes of diagrams can also be treated in collinear factorization. 
The Glasma graph would refer to double parton scattering and is suppressed by 
Aqcd/Q^ relative to a single parton scattering diagranEHlElll (such as our Jet 
Graph). The high-energy limit of QCD therefore allows us to study novel QGD 
processes which would otherwise be largely inaccessible by experiment. 

The motivation for the first study of two-particle correlations in the color glass 
condensate frameworlP^ was geared towards explaining the ridge in heavy-ion col¬ 
lisions. The basic idea is that the longitudinal chromo-electric and -magnetic fields 
created in the collision of two large nuclei radiate correlated particles approxi¬ 
mately isotropically in the azimuthal separation of the pairs. The boost invariance 
of the classical field results in boost invariant particle production naturally explain¬ 
ing the long-range rapidity correlations. Gollimation in azimuthal separation, A0 
is achieved by final-state effects such as strong radial flow due to hydrodynamic 
expansion. 

The strength of the correlated two-particle distribution can be expressed a^^ 

/ \ - ^ / dN \/ dN \ 

\ dypd‘^pj_dygd‘^(ij_ Lonn. \ dypd^pj_ / \ dyqd‘^n± / ' 

We should stress that the above is an intrinsic correlation expressed in terms of the 
single particle distributiorj^and is essentially a consequence of dimensional analysis 
- the correlations are entirely classical in nature and there is only one scale, Qs 
characterizing the strength of the fields. The factor 1/(5'±Q|) has the physical 
interpretation of the transverse area of a flux-tube from where the particles are 
emitted divided by the transverse area of the system. The coefficient hz can be 
calculated but it is a constant of 0(1) and the strength of the resulting correlation 
fares well with the experiment in a blast-wave model.l^^ The subsequent discovery 
of the double-ridge (a symmetric collimation on both the near and away-side) which 
was found after careful subtraction of the recoiling jet leaves a gap in this picture. 

It was later realized that the glasma graphs contain an intrinsic dependence on 


^In general, the double inclusive distribution can be written as 

d^N \_/ dN W dN \ / d2jv 

dypd^p±dyqd^(l± / \ dypd'^p± / \ / \ dypd^p±dyqd'^ci± 


( 30 ) 
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Fig. 30. Anatomy of a proton-proton collision. The away-side peak, associated with mini-jet pro¬ 
duction can be understood from the jet-graph (two gluons produced from a single ladder) as 
shown in the right diagram along with its schematic contribution to the per-trigger-yield plotted 
in blue. The Glasma graph contribution (left diagram) is shown schematically by the orange curve. 
The shaded gray region (extracted experimentally by the ZYAM procedure) is referred to as the 
associated yield 


the azimuthal separation irrespective of the presence of an additional collimation 
mechanisnP^ which is not apparent in equation (29). After the observation of a 
ridge in high-multiplicity proton collisions a blast-wave analysis showed that the 
near-side ridge in p-p could not be described by combination of isotropic production 
from glasma flux tubes with subsequent transverse flowP^ Instead, the systemat- 
ics of the p-p ridge were consistent with the intrinsic azimuthal dependence from 
glasma graphs leaving little room for any additional transverse expansion. 

We therefore have the following picture of the dihadron correlations in proton- 
proton collisions (summarized graphically in Fig. 30). The Jet Graph generates 
particles predominately back-to-back {i.e. at relative azimuthal angle of tt) and 
when expressed as a per-trigger-yieldis approximately independent of multiplicity, 
an observation consistent with the experimental data. 

Due to a novel interference phenomenon generated from the intrinsic correlations 
in the nucleus (the precise form of which will be discussed in the subsequent sec¬ 
tions) the glasma graph has enhanced production on angular separations of A0 = 0 
and A0 = tt. The collimation at A0 = 0 is responsible for the near-side ridge while 


•^The per-trigger-yield is the double-inclusive distribution A) divided by one power of the single 
inclusive distribution (dN). Since the jet and single-gluon graphs proceed via a single t-channel 
exchange they have the same number of connections to the larger-x sources and therefore much 
of the centrality dependence cancels in the ratio. 
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that at A0 = tt explains the double-ridge (after appropriate jet subtraction). 

The glasuia graph has a much stronger centrality dependence than the jet graph. 
In the dilute-dilute limit the per-trigger yield of the glasma graph scales with 
coupling as while in the dense-dense limit it goes as This large enhancement 
is expected for high multiplicity events at moderate ^ Qs- 

With this basic picture in mind we now go on to discuss two-gluon production in 
the color glass-condensate framework. The next two sections focus on the perturba¬ 
tive p± ^ Qs limit of di-jet production and glasma graphs respectively and explore 
the mechanism behind the intrinsic azimuthal correlations. We will then give an 
overview of the relevant phenomenology of this framework and discuss some open 
issues. 

Mini-jet production In this section we explore in more detail our understanding 
of di-jet or mini-jet production. The away-side (or recoiling jet) is pervasive in all 
dihadron analyses and serves as an important baseline. 

Computing the di-jet amplitude in the Color-Glass-Condensate framework is a 
formidable task. The diagram shown in the upper right of table involves calcu¬ 
lating the gluon propagator in the full nonlinear classical background field. In the 
dilute-dense limit this was looked at ^^d the resulting expression could 

not be written in a /^T-factorized form as the end result involved correlators of four 
Wilson lines. This quadrupole operator obeys its own evolution equation and it 
cannot be factorized into a simple product of dipoles. We point out, however, that 
its nonlinear Gaussian approximatiorP^ in terms of the dipole operator agrees 
numerically very well with the full result in JIMWLK evolution.^^^ 

In contrast, di-jet production has been extensively studied in collinear factor¬ 
ization and at leading order the partons are produced strictly back-to-bacl^ This 
is due to the fact that the internal momentum distribution of the partons in the 
hadrons is neglected and conservation of momentum requires the produced par- 
tons to have equal transverse momentum squared and opposite directions, z.e., 
a ^ ^^(pi — q^)(5(|0p — — tt). Any internal motion, for example due to Fermi 

motion or soft emissions, will lead to a decorrelation of the jets. 

The soft radiation emitted between the produced gluons and their respective 
hadrons in addition to the emission of rapidity ordered gluons between the two 


®The single inclusive spectra scale as as and a'j'^ in the dilute-dilute, dilute-dense, and dense- 
dense limits respectively. 

^The collinear factorized form for di-jet production can be obtained from the -factorized form 
by substitution of the unintegrated gluon di stri bution with the appropriate parton distributions, 
= xg{x) (k^). From equation |33| we obtain 


d?N 

dyp(Pp^dyq(P(l^ 


= Xpg{Xp) Xqg(Xq) 


4a^iV| i5^^ypx+q±) 
(JV| - 1) p2_qi 


(31) 


agreeing with equation (60) oiPMlwhen the large-r/ matrix element given in equation (54) is used. 
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Fig. 31. Two-gluon production versus relative azimuthal separation for min. bias p+p collisions 
(details are labelled in the figure). The three curves represent three approximations; Dashed blue: 
Eikonal approximation (parametrically small angle scattering of incident particles) to the 2 — 4 
matrix element without radiation between the tagged gluons (as Ay —)■ 0), Solid purple: Leading 
term in Ay —)■ oo of the QMRK result, Dotted orange; Includes the full BFKL Green’s function 
taking into account the emission of rapidity ordered gluons between the tagged final-state partons. 


tagged partons is encoded within the following kr factorized di-jet cross-section 

<PN S, 


(Pp±d^(l±dypdyg - 1) p^q5_ 


(32) 


X / $p(ko_L)$T(k3_L)f/(ko_L-p±,k3_L+q±,Ayp,). 

Jkoi,k3^ 


In the context of BFKL resummation (where the above is rigorously defined) the 
quantities ^ are the BFKL evolved impact factors and are known to next-to-leading 
order.l^SH this work we will apply a hybrid formalism that identifies the impact 
factor with the unintegrated gluon distributions introduced earlier. We employ in 
this work the LO BFKL Green’s function, the explict form of which can be found 
ii]205] example. Though its NLL form is knowrP^ the corrections introduced to 
the observables considered here are small enough to be neglected. 

The BFKL radiation can be switched off hy taking the as Ay limit of the above 
Green’s function and the resulting di-jet cross section becomes 


d^N 

dypd‘^pj_dygd‘^qj_ 


Sx_ 

87r8(iV2 - 1) plcil 


/ 


d2k_L$(k_L)i>(k_L 


+ p± + q±) 


(33) 


The above expression is valid in the multi-Regge kinematics (MRK) corresponding 
to the leading Ay oc limit of the 2^4 amplitude. Using the 2^4 matrix 
element in the eikonal approximation results in the quasi-multi-Regge kinematics 
(QMRK) whose precise form can be found in.l^SlI 
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Figure [^demonstrates the results of the above discussion and shows the relative 
importance of including rapidity ordered gluon emissions even for moderate as Ay 
as encoded in the BFKL formalism. Neglecting these emission would spoil the 
agreement with the experimentally measured away-side mini-jet. The details of 
the full calculations will be discussed further on, but The agreement between the 
BFKL evolution and the experimental data can be seen the agreement of the away- 
side yield in figure for p-p collisions. Another example is the agreement between 
the blue curve in the right plot of figure that shows the mini-jet with BFKL 
evolution in peripheral p-Pb and the corresponding ATLAS measurement of the 
same. The mini-jet itself contains interesting QCD dynamics and may be a useful 
probe of Pomeron exchange and BFKL evolutiorP^ complementing studies that 
look for the growth in di-jet cross section as proposed in.^^^ 


Glasma Graphs We now review the mechanism behind the intrinsic collimation 
that was predict ecGHHEl^ to exist in systems lacking hydrodynamic flow or final state 
rescattering. We will provide an overview of the glasma graphs that are responsible 
for this intrinsic azimuthal correlation and argue that the resulting signal is a 
quantum interference effect sensitive to correlations present in the nuclear wave- 
function before the collision occurs. 

In order to set the stage and develop some notation we will begin with an 
overview of the perturbative ^ Qs calculation of single gluon production and its 
resulting /cx-factorized expression. The starting point is the probability to produce 
a gluon having rapidity y and transverse momentum k^, 

^ .^(|Mnk)P). (34) 


dyd‘^kj_ 2{2 'kY 


a, A 


The classical contribution to the amplitude for the production of a single gluon 
with four-momentum k is Mx(k) = ej^(k) k‘^A^(k) with representing the gluon 
polarization vector having polarization states A. The above should be evaluated in 
the on-shell k^ ^ 0 limit. The expectation value (• • •) represents an averaging over 
the color configurations of the two nuclei. 

At large transverse momentum the gauge fields can be expressed explicitly in 
terms of color charge densitie^lHllll221 of the projectile and target nuclei, pp and pp , 




-fabCg3 


/ 


(27r)‘ 




Pp4k±)Pr(P-L -k±) 

ki(p_L - k_L)2 


(35) 


L^{p±, kj_) is the effective Lipatov gluon emission vertex?^ The resulting expecta- 
tion value in equation (34) is a product of two-point correlators over the color charge 
densities of the target and projectile. In the non-local MV model the momentum 
space correlator is 

(/3“(kii)p*4k2±)) = (27r)2<i“V"(ki±)<^(') (kix - k2^) ■ (36) 

We should stress that in order to accommodate the small-x evolution one needs to 
generalize the local MV model of equation (23) to include a nonlocal version such 
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that — y±) and consequently in momentum-space a function of the 

transverse momentum (z.e. /i^(|kx|)). Further discussion of this subtlety can 

be found inP^ Combing equations (28), (34), (35), (36) the following /^T-factorized 


form for the single inclusive gluon production is obtained 
dN asNc Si_ f 


/ 


$p(kx)$,(p±-k^). 


(37) 


dypd'^pj_ - 1) p\ J (27r)2 

We should point out that the above factorized form is obtained in a different limit 
than the usnal-kr factorization derived from the scattering of a dilute projectile off 
a quasi-classical target as examined irP^ for example. Equation ^ has been able 
to successfully describe a wealth of single-inclusive data from both p-p and d-Au 
collisions \2121\21M 

We now consider the production of a pair of gluons with transverse momenta 
Pj_ and and longitudinal rapidity i/p and i/q within the same framework. The 
double-inclusive gluon distribution is 

d^N 1 


dypd?p^dyqd?(\^ 4(27r)^ 


E (I^A'(p,q)P) • 


(38) 


a,b,X,a 


The classical contribution to the two-gluon matrix element is M^^(p,q) = 
ej^(p)e^(q) p^A^’^(p) g^A^’^(q) and these gauge fields can be expanded in terms 


of the color charge densities when px,qx ^ Qs as done in equation (35). Putting 


all the pieces together we obtain the following schematic form of the inclusive two- 
gluon distribution 

d^N 

d^pd^pxd^gd^qx 


[ _7'(p± , q± , kjX ) (Pp (ki_L ) Pp* (kgx )Pp (k2_L )/?; (k4X ) ) 


X (Pt(P-L - kix)Pp (P-L - k3x)/?p(qx - k 2 x)Pp (q-L - k4x)> 


where color indices have been suppressed and a shorthand notation, T, has been 
used for the product of multiple Lipatov vertices, the detailed form of which will 
not be needed for the qualitative discussion to follow. 

The basis of the glasma graph framework is the factorization of the four-point 
functions into products of two point functions using the Gaussian form of the color 
charge distribution such as the local MV model in equation ( [^ or the non-local MV 
model in equation ( [^ (used for the remainder of this section) that accommodates 
small-x evolution. The resulting factorization results in a total of nine diagrams, 
each carrying a unique dependence on the azimuthal separation of the two produced 
gluons which we will discuss in turn. 

An important piece of the above formalism, which remains to be discussed, 
is the rapidity dependence of the above correlation. At leading order in ag {i.e. 
working within the local MV model) the particle production is boost-invariant and 
therefore independent of the rapidity separation of the pairs and their respective 
separation from the leading hadrons. Quantum corrections modify this picture and 
lead to the breaking of boost invariance. The difficulty in the proper treatment of 






















50 Kevin Dusling, Wei Li, Bjorn Schenke 


the above problem stems from the fact that not only must quantum evolution be 
included between the measured gluons and the leading nuclei but also between the 
two measured gluons. A dense-dense factorization formalism resumming leading 
logarithms in 1/x corrections for inclusive observables developed irPI^H^I^was used 
to compute the rapidity dependence of gluon-pair production in.^SHl However, this 
factorization may break down when ^ Regardless, the framework used 

here serves as a useful starting point for addressing the role of long-range rapidity 
correlations in the dense-dense limit. 

As noted above the glasma graphs are obtained by forming all possible contrac¬ 
tions over the color sources. One of the nine resulting graphs corresponds to the 
single-inclusive spectra squared. While it is leading in Nc it has no dependence in az¬ 
imuthal separation within the approximations used here. However, event-by-event 
fluctuations or an impact parameter profile of the nuclei would indeed generate 
azimuthal correlations within this framework (as seen for example irP^. We now 
focus on the remaining eight diagrams and the azimuthal correlations they generate 
without introducing a global geometry. 

Let us give as an example the expression from one such diagram responsible for 
the near-side ridge 


d'^p±d‘^ci±dypdyg 


f 

47ri0C {N! - 1)3 piqi L 


$p(kx)$,(px 


k±)$T(q-L-k_L) , 


where the unintegrated gluon distributions are evaluated at (small) x values on the 
order x ^ pt j\fse^^^ - the precise prescription is given inESIlISSSl The feynman 
diagram is shown in figure There is a single loop momentum in this diagram 

that is integrated over. The transverse momentum flowing out of the nuclear 
projectile and target (represented here by the blobs) is of the order of their respec¬ 
tive saturation momentum. In this case we therefore must have \p±_ — k^| ^ Qs 
and \q_j_ — k^l ^ Qs- Furthermore the magnitude of the loop momentum must 
also be near the saturation scale |kx| ^ Qs- This leads to the conditions that 
q]_ — 2q±Qs cos (l)kq = 0 and p^ W PtQs cos (l)kp = 0. This pair of equations con¬ 
strains the produced gluons such that 1) (^p = collimation on the near side. 2) 
p±_ ^ maximum correlation for similar transverse momentum and 3) p±_ Qs 
maximum correlation near a semi-hard scale. All three of these features are born 
out by the full numerical calculation and consistent with the trends seen in the 
data. 

The full glasma graph calculation includes contributions from all eight graphs - 
the detailed expression can be found in.All eight diagrams provide a near- and/or 
away-side collimation and are necessary for a quantitative description of the data. 
Two diagrams of note are those having delta function correlations, d?N/d/S.(j) ^ 
± q^), which was interpreted as HBT-like correlations in.l^^In the dilute- 
dense limit (by dilute we mean (/)p(kx) ^ iFese two graphs are the only 

glasma graphs that remain and are therefore crucial for our understanding of highly 
asymmetric collisions. The delta functions will be broadened by hadronization. For 
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Fig. 32. Momentum flow demonstrating the generation of the near-side collimation. The intrinsic 
parton momentum of either hadron is on the order Qs and therefore |p_L — k_i_ | ~ |q_L — k_L | ~ Qs 
where k^ is a loop momentum constrained such that |k^| ^ Qs- 


practical considerations the delta function is replaced by a Gaussian, the precise 
form of which can be found in,^^^ originally fit to the p-p dihadron correlations 
before the availability of any data on p-Pb. The good agreement with the subsequent 
p-Pb data shows the robustness of this modeling of non-perturbative hadronization 
dynamics. We should stress that the associated yield (the integral over on the 
near side) is insensitive to the functional form of the smearing function. 

In order to convert the two-gluon correlations presented above into the hadronic 
observables, an appropriate hadronization procedure must be used. In what fol¬ 
lows we show results using the NLO KKP parameterizatiorP^ of fragmentation 
functions for gluon to charged hadrons. It has recently been found that the NLO 
KKP results are troublesome at LHC energie^^^ and suggested that the gluon-to- 
hadron fragmentation functions were a probable source of this problem. Extraction 
of new fragmentation functions by fits to the more recently available single inclusive 
hadronic data would be highly valuable. 

Figure shows a comparison of the di-hadron correlation compared to the 
results for high multiplicity p-p collisions. The full numerical calculations of all eight 
glasma graphs along with the away-side mini-jet contribution with BFKL evolution 
corroborate the qualitative picture shown in figure The centrality dependence 
is controlled by an appropriate choice of initial saturation scale Q^q that fixes the 
initial condition in the rcBK evolution equation. Fits to deep-inelastic scattering 
constrain = 0.168 GeV^ and we take this value as representative of min. bias 
p-p collisions. For convenience, we work with integer multiples of this saturation 
scale. For example, central p-p > 110) corresponds to 5-6 times this min. 

bias value. 
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Fig. 33. Long range (2 < A 77 < 4) per-trigger yield of charged hadrons as a function of Acf) for 
p-p collisions at \/s = 7 TeV. Data points are from the CMS collaboration. The curves show the 
results for Qq{x = 10 “^) = 0.840 GeV^ and Qq{x = 10 “^) = 1.008 GeV^. 



Fig. 34. Representative sample of p-Pb data form the CMS, ALICE and ATLAS collaborations of 
the per-trigger yield along with calculations within the glasma graph framework. 


The overall strength of the glasma graph contribution is controlled by ag which 
is evaluated at one-loop running at the relevant momentum scale of the process. 
In addition, a correction from non-perturbative dynamics (for example multiple 
scattering) is taken into account through the multiplicative pre-factor l/( that 
enhances the glasma graph contribution relative to the jet contribution. 

The parameter ( is particularly sensitive to the multiplicity distribution (see 
figure and independent fits corroborate the value of C = 1/6 used in the glasma 
graph ridge analysis. Lattice calculation^li^UllI find that this constant can be small, 
C ^ 0.2 — 1, lending support that non-perturbative corrections due to multiple¬ 
scattering enhance the signal. 

A comprehensive comparison of the glasma graph framework with all the avail¬ 
able p-p, p-A and d-A data was presented irP^ and will not be reproduced here. 
Instead in figurewe show a representative set of p-Pb data from the CMS, ALICE 
and ATLAS collaborations along with the corresponding glasma graph calculations. 

In summary, the glasma graph framework is able to account for many features 
of the data on a qualitative and quantitative level. These include 1) the long range 
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nature of the correlations, 2 ) the nearly symmetric near- and away-side ridge (the 
double ridge), 3) the strength of the correlation as a function of transverse momen¬ 
tum, in particular the maximum correlation for px ^ Qi. ^ Qs? 4) the non-trivial 
centrality dependence of both p-p and p-Pb collisions. 

Since the time of the above analysis more precise data on the experimental side 
have uncovered two shortcomings in the glasma graph description. The first is the 
presence of a large and positive in central p-Pb collisions. The glasma graphs only 
contribute to the even harmonics and the peripheral jet diagram gives a negative 
V 2 ,. Whether modifications to the jet can accommodate a positive has yet to be 
fully explored. Further, the interference diagram between the jet and glasma graph 
also appears to vanish.l^^ Recent calculations within classical Yang-Mills provide 
one source of positive (see section 4.4) but this has yet to confront the data. 

A final shortcoming of the glasma graph description (at least within the frame¬ 
work used here) is the negative C 2 { 4 } measured by CMS. The glasma graph and jet 
graph both give a positive contribution to C 2 { 4 }. The empirical observation that 

—C 2 { 8 } suggests that the correlation is driven by corre¬ 
lations among the disconnected diagrams at the single-inclusive level. We should 
stress that the multiple particle correlations are dominated by soft momentum 
P_L ^ Qs where the glasma graph framework no longer applies. It would be very 


interesting to see if and where C 2 { 4 } changes sign as the minimum pt used in the 
analysis is increased and when a rapidity gap between the particles is introduced 
in order to suppress jet-like correlations. 


Multi-gluon production in the CGC The framework laid out for two-gluon 
production in the previous section has been extended to three gluonJ^^ and gen¬ 
eralized to n-gluons.l^^ 

The results can be represented diagrammatically as glasma graphs in a similar 
fashion as above and the subset of these diagrams that contribute for p± ^ Qs- 
has been identified and the n-gluon distribution is 

/ d-N \ ^ C(A^,^-1) Q%S^ / T 

\dyid‘^kij_---dynd‘^kn±/ ' 27r ■ ■ ■ k^_L 

(39) 

where the sensitivity to the infrared dynamics is encoded in the constant ( intro¬ 
duced earlier. The factorial cumulants (d^N) are those of a negative binomial 


pr(n,fc) 


r(/c + n) 
T{k)r{n^l) 


(40) 


where n is the average multiplicity and k = C{^c ~ l)Q|'S'x/(27r) is proportional 
to the number of flux tubes. For k = 1 the negative binomial is a Bose-Einstein 
distribution and for k ^ oo the distribution is Poissonian. The left plot in figure 35 
shows the resulting multiplicity distribution in p-p collisions in the IP-Sat modeP^ 
with the one free parameter ( = 0.155 fit to data at 0.9 TeV and then exhibited for 
other energies. 
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Fig. 35. Multiplicity distribution within A:_l- factorization in the IP-Sat mod el^^ (left) and from 
classical Yang-Mills with intrinsic fluctuations of the proton scaturation scal d^^^ l (right). 


The right plot in figure 35 shows the recent results oP^ for the multiplicity 
distribution in p-p from classical Yang-Mills simulations. In this case there is no 
free parameter ( and the infrared behavior is instead determined by the resulting 
non-linear dynamics of the gluon fields. In this framework there are multiple sources 
of fluctuation. First there is a geometric dependence from the saturation scale in 
the IP-Sat model which sets the initial field configurations. Second, sub-nucleonic 
fluctuations from the sampling of different configurations of color charges according 
to the MV model give rise to the negative binomial distribution of gluon number 
fluctuations for fixed geometry. 

However, convolution of these two sources alone are not sufficient to account for 
the broad width of the multiplicity distribution. In addition there must be intrinsic 
fluctuations of Qs in the proton. Modeling the distribution of saturation scales of 
the proton as a Gaussian distribution provides a good description of the data. The 
key point for this discussion is that high multiplicity events easily accommodate 
saturation scales 5-6 times the min. bias values as used to explain the ridge data. 


4.3.2. Classical Yang-Mills dynamics 

The calculations within the glasma graph framework presented in the previous sec¬ 
tion are strictly valid for Pt ^ Qs and are infrared divergent if extended to lower 
transverse momentum. The divergences are regulated by the non-perturbative dy¬ 
namics of the gluon fields at scales of order Qs. These dynamics can be captured 
by solving the classical Yang-Mills equations numerically in the forward light-cone 
and are valid at any transverse momentum - only limited in the ultraviolet by the 
inverse lattice spacing and in the infrared by the inverse nuclear size. The pertur¬ 
bative, high pt, limit of the classical Yang-Mills simulations can be expressed in 
terms of similar glasma graphs if a Gaussian distribution in the sources is used. Gor- 
rections due to small-x evolution in numerical simulations of quasi-classical gluon 
production have been dond^^ but in the following section on classical Yang-Mills 
dynamics a local Gaussian distribution of sources is used (which does not accom- 
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modate small-x evolution). 

Empirically, the near-side ridge persists to momenta well below the saturation 
scale and it is therefore imperative to investigate whether the Glasma flux-tube 
picture holds in this kinematic region. It is also useful to establish the size of 
the corrections due to multiple scattering at larger transverse momentum. In the 
remainder of this section we will review the numerical result^^^ of the classical 
Yang-Mills simulations. We will discuss how the classical Yang-Mills results com¬ 
plement the glasma graph description and also provide possible resolutions to some 
of the discrepancies with data alluded to in the previous section on glasma graph 
phenomenology. 

In the CGC picture, the dynamics of a high energy collision is described by the 
solution of the classical Yang-Mills equations 

(41) 

in the forward light-cone of the collision. In the above expression is an eikonal 
color current generated by a color neutral distribution of classical charges moving 
along the light-cone, 

= S^+pp{x~,x±) + 6^~pj.{x''~,x±) . (42) 

The sources are distributed according to a Gaussican distribution having variance 

5^ {p°'{x±)p\y±)) = 5'p/Pb(b_L)V'’(5(2)(x_L - y±) . (43) 

The function 5'p/pb(bx) encodes information on the spatial structure of color charge 
distributions of the projectile and target. The results shown here use a constituent 
quark proton model whereby the color charge density is concentrated around the 
transverse positions of three constituent quarks whose positions fluctuate event to 
event. The use of a spherical proton has a negligible effect on the results.l^^ The 
result is even less sensitive to the impact parameter dependence of the color charge 
in the nucleus which is sampled from a Wood-Saxon distribution. 

Starting from field configurations derived from the color charges in a given event 
the classical equations of motion are solved and the gluon distribution extracted 
at a given proper-time by measuring equal-time correlation functions of the gauge 
fields. Using either two-particle correlations or the event plane method V 2 {pt) and 
V 3 {pt) can be extracted from the gluon spectrum as show in the left and right plot 
of figure We now discuss the implications of these results. 

It is striking that a sizable V 2 {pt) is present at the initial time r = 0+ as it 
therefore cannot be attributed to a collective expansion of the system. This is in 
contrast to a hydrodynamic interpretation where the elliptic flow is initially zero and 
gradually develops over the lifetime of the collision. Further evolution of the Yang- 
Mills simulation finds only a modest decrease in the V 2 at later times indicating 
that the glasma flux-tube picture is robust against gluon re-scattering (a similar 
conclusion was reached by earlier lattice simulationJ^^ as well). 
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Fig. 36. (Color online) Gluon V2(pt) (left) and V 3 (pt) (right) at times r = 0, 0.2, 0.4 fm/c in p-Pb 
collisions at impact parameter = 0 in the constituent quark proton model. Open and closed 
symbols correspond to results obtained using the event plane method and two particle correlations, 
respectively. Error bands include statistical errors only. Experimental results by the ATLA^^and 
CMS collaboratior®! for inclusive hadrons are also shown as a guideline for comparison. 


The observation that V 2 as extracted from the event-plane method and two 
particle correlations are similar within errors suggests that the dominant origin of 
the V 2 is the breaking of rotational symmetry of the single particle spectrum on a 
per-event basis. To elaborate on this point, the initial-state {i.e. r = 0+) classical 
Yang-Mills results can be divided into two contributions. One consisting of glasma 
graphs corresponding to the connected diagrams shown in the previous section. The 
second contribution stems from seemingly disconnected graphs (disconnected on a 
perturbative level) but connected when the two gluons scatter on a common color 
field domain. This is similar in spirit to the effect seen iipMl (where the proton was 
treated as a dilute object) and will be discussed in further detail in the upcoming 
section. 

The classical simulations fail to reproduce the magnitude of the V 2 {pt) at larger 
values of transverse momentum. The situation should worsen when hadronization 
is taken into account. Hadrons at a given px fragment from gluons of higher px and 
therefore the V 2 shown in figure [36] will be rescaled to lower px when hadronization 
is accounted for. This demonstrates the importance of including small-x evolution 
in the nuclear wavefunction absent in the classical Yang-Mills simulations. The clas¬ 
sical Yang-Mills and glasma graph picture complement each-other quite nicely; the 
latter responsible for the large Pt ^ Qs ridge due to intrinsic {i.e non-factorizable) 
two-parton correlations including small-x evolution and the prior responsible for 
the softer ridge, Pt ^ Qs generated from the event-by-event breaking of rotational 
invariance. 

We now turn to a discussion of the gluon 'r’ 3 (px) as shown in the right plot 
of figure 36 The first observation is that the initial 'T’s at r = 0+ is zero. This 
is a consequence of the initial gluon spectra being symmetric under kx ^ — kx 
resulting in the odd Fouier harmonics being identically zero at the initial time. This 
is consistent with the vanishing of odd harmonics in the glasma graph calculation. 
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At a proper-time ofr = 0 . 2 fm/ca sizable has built up. Further evolution to r = 
0.4 fm/c, where the system is practically free-streaming, results in a modest increase 
in the signal. The agreement between the event plane method and two particle 
correlations suggests that the event-by-event breaking of translational invariance 
combined with coherent final state effects is ultimately responsible for the v^. The 
most surprising finding is the lack of correlation between the global initial state 
eccentricity, 63 , and the px integrated on a per-event basis; refuting the notion 
that the source of ^3 is from the build up of a global energy-momentum flow via 
classical Yang-Mills dynamics. 

Finally, we stress that the Yang Mills calculatiorP^ is not able to reproduce 
the azimuthal anisotropy coefficients in A-A collisions, because of the large amount 
of uncorrelated color field domains present in a single collision and the lack of 
hydrodynamic evolution. This further confirms the standard interpretation of the 
Vn in heavy ion collisions as strong final state effects. 

To summarize, the azimuthal correlations generated via the classical Yang-Mills 
simulations can be separated into three sources: 1 ) A genuine non-factorizable two- 
particle correlation - the glasma graph contribution 2 ) perturbatively disconnected 
graphs connected by the event-by-event breaking of rotational symmetry and 3) A 
contribution from final state interactions generated by the Yang-Mills dynamics. In 
the following section we review work that has focused on the second contribution; 
the event-by-event breaking of rotational symmetry. 


4.4. Scattering from color-field domains 

The basic idea can be traced back to Kovner and Lublinsk 3 EMl[SSll who argued 
that domains of color-electric fields of size l/Qs will produce non-trivial angular 
correlations in the single-particle distribution. Following their arguments, this can 
be easily seen by considering the scattering of two projectiles from a stationary (in 
the lab frame) hadronic target. The two projectiles may have different rapidities 
but at high enough energies their wave functions are boost invariant. Therefore 
the long-range nature of the correlations are a general consequence of the high- 
energy limit. If the two projectiles strike the target at the same impact parameter 
both will see the same color configuration in the target and be scattered in similar 
directions. This provides the source of near-side correlations. If the projectiles under 
consideration are gluons, it is just as likely that they can have opposite charge (as 
gluons transform in the adjoint representation of SU(3) ) and therefore scatter 
in opposite directions. For quarks in the fundamental representation this would 
not be the case and the near-/away-side symmetry would be broken. As the size 
of correlated domains in the target are of order l/( 5 |(target) we parametrically 
expect the strength of the two particle correlation to go as (Q|5'x) where 5'x is 
the overlap area. 

So far this is conceptually identical to what has been done in the previous section 
on glasma graphs except considered in a different Lorentz frame. Indeed the work 
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Fig. 37. Left: Schematic illustration of localized domains of chromo-electric fields in the tar¬ 
get. Right: The dipole operator 1/NcTr R(x^)yt'(y^) after evolutio n by the JIMWLK equation 
demonstrating the persistence of local fiuctuations of the Wilson line.E^ 


opm obtains the same glasma graphs when using a Gaussian weight functional for 
the color charge densities. 

It was argued however that the Gaussian averaging procedure may underesti¬ 
mate the strength of the correlation. In the Gaussian approximation the two-gluon 
production is written as a term that factorizes into a product of single inclusive dis¬ 
tributions and eight remaining connected diagrams (the glasma graphs considered 
earlier), the latter of which are suppressed, 

/ \ / dN \ / dN \ 1 ^ 

\ dypd^p±dygd^(ii_ )\ dypd?p^ ) \ dygd^q ^® ® 

(44) 

However, in the Gaussian approximation this first term has no azimuthal correla¬ 
tions - it corresponds to the independent emission of two gluons. An analysis going 
beyond leading order in Nc has found that the four-point function does not nec¬ 
essarily factorize into a product of two-point function^^^ and these factorization 
breaking terms contribute at leading order in Nc to correlated two-gluon produc¬ 
tion. 

A model for the generation of initial state azimuthal asymmetry from scattering 
from color domains was introduced inl^^It was shown that scattering of a dipole of 
size rsj I Ipt from a target with fixed configuration of color electric field reproduces 
the pt dependence of V 2 and V 4 for a sufficiently polarized target. 

Gonventionally in the MV model the averaging over the field configurations of 
the target using a Guassian correlator, 

£ {EnW)E^{h2)) = (jv/- - b2) (45) 

will result in isotropic particle production. In order to take into account observables 
that are sensitive to the breaking of translation invariance on a per-event level a 
model is introduced such that the target configurations are constrained such that 
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Fig. 38. Comparison of the color domain model oB^with r’ 2 (p±) (left) and r’ 3 (p_L) (right) with 
high multiplicity p+Pb collisions. 



the electric fields point in a specific direction d, 

£ {Et{h,)E^,{h,)) = -b2) 


“I" ^“4 


CLiCij 



(46) 


The above expression takes into account that in the vicinity of bi ^ b 2 there is 
a subclass of events that have a predefined orientation and the average over these 
subclasses are done after the observable has been computed. 

In this framework there exists a correlation even at the single particle level, 


(cos n{(j)p - (/)a )) (p±) = J d(l)p (cos n{(pp - (ps. ) 

and will correspondingly produce azimuthal correction in the disconnected terms 
in higher order cumulants. 

One of the strengths of the above model is that it is able to accommodate a 
negative C2{4p^ 






4(7V|-1)3; 

Furthermore, the consideration of C-odd fluctations in the target, related to the 
odderon, are able to generate a V 3 . 

While the above color domain model appears to use only Gaussian fluctuations, 
(i.e. the only correlator used equation (46) is a two-point function) the separation 
of averages is outside of the assumptions of the CGC framework. Implicit in this 
separation of averages is that the angular fluctuations of the large-x color sources, 
rdo, evolve on a slower times scale. First the fast modes are averaged over with the 
modified two point function equation (46), second the observable is computed, and 
only then is the averaging over all possible electric field directions, d performed as 
required by Gauge invariance. 

A recent worlP^ has made the connection between the color domain model 
and non-gaussianites in the color glass condensate. We direct the reader there for 
further discussion of the interpretation of the polarization factor A in the color 
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domain model and its relation to the possible sources of non-Gaussianity in the 
CGC effective theory. 

4.5. Conclusions on initial state interactions in small systems 

In conclusion, the previous sections have demonstrated how glasma graphs and 
classical Yang-Mills complement each other as two different approximations of the 
same effective theory of high energy nuclear collisions able to address different fea¬ 
tures of the data. The connection to the scattering from localized domains of color 
electric fields with finite polarization can be made by considering non-Gaussianities 
of the color charge correlators. 

We conclude this section by briefly discussing two other initial state proposals. 
While these have yet to address the data at the same quantitate level as the hydro- 
dynamic and GGC models discussed thus far, these proposals contain interesting 
ideas which should be explored further. 

The work demonstrates that angular correlations can be generated by 

the exchange of two BFKL Pomeron ladders. The transverse momentum exchanged 
by the Pomeron is conjugate to the impact parameter of the collision and as events 
with a fixed multiplicity correspond to a finite impact parameter, an asymmetry is 
generated. Motivated by Gribov Reggeon and Pomeron calculus it was shown that 
the correlation is long-range in rapidity and that a double-ridge signal is generated. 

A more recent worlP^ argued that non-abelian beam jet bremsstrahlung gen¬ 
erates long range rapidity correlations having a hierarchy of non-trivial azimuthal 
harmonics. The anisotropies are generated from radiating clusters that have accu¬ 
mulated a net transverse momentum kick on an event-by-event basis. In addition 
to explaining the ridge in small colliding systems, it was argued that this mech¬ 
anism may be able to explain the approximate energy independence of azimuthal 
asymmetries observed in the RHIC beam energy scan. 


5. Future Directions 

While a lot of progress has been made towards understanding novel long-range cor¬ 
relation phenomena in small colliding systems, many questions remain open, requir¬ 
ing future efforts by both experimental and theoretical communities. We consider 
the following two big questions we feel should be addressed in order for progress 
to be made in the coming years: 1) Have we reached a consensus on the nature of 
long-range correlation phenomena in small (and/or large) systems? 2) If a strongly- 
coupled QGP is indeed formed, what fundamentally new knowledge do we gain from 
small systems? 

To answer these questions, we discuss a few possible future directions in this 
section: 

• Further scrutiny of hydrodynamic paradigm 

• Nature of the ridge in p-p collisions 
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Fig. 39. Cumulant V 2 ratios, V 2 {^}/v 2 {^} (left) and r’ 2 { 8 }/'^ 2 { 6 } (right), as a function of 
'^ 2 { 4 }/r’ 2 { 2 } ratios for high-multiplicity events in p-Pb collisions at = 5.02 TeV and Pb- 

Pb collisions at = 2.76 TeV.^^ The solid curves are predictions from a proposed universal 

behavior of fluctuation-driven eccentricities.^^ 


• Jet-medium interactions in p-p and p-Pb collisions 

• Pre-equilibrium dynamics 


5.1. Further scrutiny of hydrodynamic paradigm 

The scrutiny of standard paradigms by identifying inconsistencies with experimen¬ 
tal observations has been the recipe for many major scientific discoveries. In this 
light, precision studies of the hydrodynamic framework in small systems should be 
steadily pursued. The main theoretical challenge for p-p and p-A systems lies in 
the large uncertainties in determining the initial-state geometry, due to the lack of 
knowledge of the proton structure as well as its fluctuations, as discussed in Sec. 4.1. 
To make progress, identifying observables that are primarily sensitive to either the 
initial-state geometry or final-state dynamics is the most promising approach. 

It was argued in Ref.,1^ that for purely fluctuation driven initial-state 
anisotropies, such as in p-p and p-A collisions, the initial-state eccentricity distri¬ 
bution is universal. With the assumption of V 2 ^ 62 motivated by hydrodynamics, 
the V 2 cumulants measured from 2-, 4-, 6 -, and 8 -particle correlations are predicted 
to follow a specific relation, as shown by the solid lines in Fig. for ratios of 
'^ 2 { 6 }/'i; 2 { 4 } and V 2 {^}/v 2 { 6 } as a function of V 2 {^}/v 2 {‘ 2 }. The recent experimen¬ 
tal results in p-Pb collisionJ^ (also shown in Fig. [3^, seem to favor the theoretical 
predictions but clearly improvement in the experimental precision is needed before 
firm conclusions can be drawn. Additional information on the initial state fluctu¬ 
ations can also be extracted by measuring event-by-event V 2 distribution in p-Pb 
collisions, as was done in Pb-Pb collisions 

Correlations among flow harmonics and the event planes (T^) are another 
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class of observable that may provide a valuable way of disentangling the effects 
of initial-state geometry from final-state evolution within hydro dynamical models. 
Such correlations were studied in great detail by the ATLAS collaboration in Pb-Pb 
collisions Emmi 

Similar measurements should be pursued in p-Pb collisions and compared to 
hydrodynamic calculations. For example, the correlation between V 2 and U 4 , or T 2 
and T 4 , is generated by the non-linear response of hydrodynamic evolution, and 
one expects such a correlation to persist in small systems as well. The interpretation 
is different for T 2 and T 3 . Final-state hydrodynamic evolution does not induce a 
strong correlation between T 2 and T 3 and indeed this correlation was observed to be 
small in Pb-Pb collisionsHowever, it has been argued that a sizable correlation 
between T 2 and T 3 may be induced by large initial-state fluctuations in small 
colliding systems such as p-Pb.^^^ Similar event plane correlations measurements 
in p-Pb could therefore provide valuable information on the initial-state geometry. 

Another observable that is mainly sensitive to the details of initial-state fluctu¬ 
ations is the transverse momentum dependent event plane fluctuations (of the same 
order). The measurement was carried out by the CMS collaboration in p-Pb and 
Pb-Pb collisions .1^ Hydrodynamic calculations have indicated that this observable 
is sensitive to the granularity of the initial energy density fluctuations, instead of 
the 77 /s of the medium .1^ Evidence of significant rapidity-dependent event plane 
fluctuations in p-Pb and Pb-Pb collisions is also observed,!^ although it may be 
sensitive to both initial- and final-state effects in the longitudinal direct ion.De¬ 
tailed studies of anisotropic flow observables and their correlations among each 
other in small systems, could provide crucial information on the initial-state in 
small systems if hydrodynamics is the correct framework to model these systems. 


5.2. Collectivity in p-p collisions 

While there has been a flood of new results on the long-range correlations in p- 
Pb collisions as summarized in Sec. little progress has been made in addressing 
the nature of the ridge in p-p collisions since its first observation in 2010. Only 
recently have new measurements of the ridge in p-p become availabld^^ with the¬ 
oretical analysis closely followingA more recent work has attempted to extract 
the V 2 harmonics of long-range correlations by a careful fitting procedure between 
peripherial and central event classes in 13TeV and 2.76 TeV.^ 

As discussed in Sec. |4.H hydrodynamic descriptions of small colliding systems 
are highly limited by the lack of knowledge of the initial conditions. Especially 
for p-A collisions, the initial-state eccentricity is highly sensitive to the shape of 
the proton and its fluctuations on very short timescales. Precision measurements 
of Vn harmonics from long-range correlations in p-p collisions should provide new 
constraints on the proton shape, and thus have the promise of improving the hy¬ 
drodynamic modeling of p-p and p-A systems. Even more precise and independent 
measurements of the proton shape and its fluctuations would be possible with an 
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Fig. 40. Left: The 1-D long-range azimuthal (A0) correlation functions in high-multiplicity p-p 
collisions at = 7TeV and p-Pb collisions at = 5.02 TeV.^ Right: long-range near-side 

associated yields for 1 GeV/c < pT < 2 GeV/c as a function of multiplicity in p-p, p-Pb and Pb-Pb 
collisions 


electron ion Collider 

The A0 correlation for hadron pairs with IGeV/c < px < 3GeV/c in p-p and 
p-Pb collisions are directly compared in Fig. (left) at a comparable multiplicity. 
While the magnitudes of the away sides are similar, the near-side ridge yield in p-p 
is about four times smaller than in p-Pb. Extracting V 2 from p-p is difficult because 
of the larger relative contribution from the away-side jet yield. The near-side ridge 
yield as a function of multiplicity in p-p, p-Pb and Pb-Pb collisions is shown in 
Fig. [4Q| (right). Above around Wrk ^ 40, the ridge yield increases roughly linearly 
with multiplicity for all systems. At a given track multiplicity, the ridge yield in p-p 
collisions is roughly 25% and 10% of those observed in Pb-Pb and p-Pb collisions, 
respectively. 

Associated yields that are collective are expected to grow linearly with event 
multiplicity, while short-range few-body correlations would be more or less indepen¬ 
dent of multiplicity (except for the bias towards enhanced jet correlations imposed 
by the selection of high multiplicity. However, this bias grows much more slowly 
than linearly with multiplicity). Therefore, pushing to sufficiently high multiplicities 
(e.g., Wrk > 150-170), the collective component of correlations would eventually 
dominate. The application of the jet-yield subtraction procedure will be more reli¬ 
able in that regime. 

To shed further light on the situation in p-p collisions, one should aim to mea¬ 
sure V 2 with multi-particle correlations. V 2 from multi-particle cumulants is less 
susceptible to jet correlations and could provide further insights. It is possible that 
a non-negligible amount of jet correlations may still be present in four-particle 
correlations for very-high-multiplicity p-p events. In that case, new methods of im¬ 
plementing an T] gap among four particles is worth pursuing to suppress short-range 
jet correlations. 

At the top LHG energy of ^/s = 13-14 TeV for p-p collisions, a data sample with 
an integrated luminosity of 50 pb“^ should be large enough to achieve the goals 


















64 Kevin Dusling, Wei Li, Bjorn Schenke 



—I—I—I—I—l— I—I— I—I—l—I— 

•°.o - 

• o 

• o 

■ o 

■ o . 

m- o 

•- -o- - 

1 1 1 1 I 1 1 1 1 I 1' 

;® NS‘2 220 


1 < < 3 GeV, 2 < IAtiI < 5 

- 

ATLAS p+Pb 

V2 

' I I I I I 

Vs^ = 5.02 TeV ; 

28 nb'^ 

. . . , 1 , ■ 

D 5 

10 

[GeV] 


Fig. 41. Left: The 1-D long-range azimuthal (A0) correlation functions for high-pT (9-12 GeV/c) 
trigger particles in pPb collisions at = 5.02 TeV. Right: The V 2 values from long-range 

correlations as a function of pT in high-multiplicity pPb collisions, up to very-high-pT region, 
with (closed) and without (open) correcting for back-to-back jet correlations.^^ 


outlined above. Cross sections of high-multiplicity p-p events increase by more than 
a factor of 10 from ^/s = 7 TeV to ^/s = 13 TeV. However, it is necessary for these 
data to be delivered under a low pile-up condition (e.g., an average pile-up of 1-2 at 
CMS and ATLAS) so that the experiments are able to trigger on high-multiplicity 
events from single p-p interactions. A special run like this lasting for a few days at 
the LHC in the future would yield very exciting physics. Novel ideas of triggering 
on high-multiplicity events under high pile-up conditions is another opportunity 
where progress can be made. 


5.3. Jet quenching in small systems 

If the observed ridge-like correlations are related to strong final-state rescatterings 
inside the medium, interactions between high-px partons and the medium should 
also be present, leading to the “jet-quenching” phenomenon first observed in heavy- 
ion collisions at RHIC. 

Due to a much smaller system size, the average path length of a parton traversing 
through a p-p or p-A system is much shorter. For this reason one might argue 
that little jet quenching should be expected. However, at similar multiplicities, a 
smaller system represents a higher energy density or temperature state. The parton 
energy loss depends on both the path length, L, and the transport coefficient 
the mean transverse momentum squared accumulated by a hard parton per unit 
length. While the average L is reduced in small systems, q increases with T^. A jet 
might lose just as much energy in a smaller but denser p-Pb system as in a larger 
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(a) PbPbv/^=2.76TeV 
= 150^b■^ 


CMS Preliminary 
0-30% centrality 


^ 26 


> 20 GeV/c 


= < 2 GeV/c 


“soft” VS “hard” 



Fig. 42. Left: The 2-D Arj-Acf) correlation function for high-pT (> 20GeV/c) trigger particles in 
Pb-Pb collisions at = 2.76 TeV. Right: The V 2 valu es at high pT 15GeV/c) versus low 

prp (rvj 1 GeV/c) for different centralities in Pb-Pb collisions .1^^ 


but more dilute Pb-Pb system. While early model calculations predicted a large 
jet quenchin^^^ one would like to see modeling with the more realistic space-time 
evolution predicted by hydrodynamic calculations. 

Initial results of high-px jets or hadrons in p-Pb collisions at the LHC indeed 
suggest no significant modification of high-px parton prodnctionESlISISl However, 
what has been explored so far mainly concerns well-reconstructed jets at very high- 
Px (e.g., > 120 GeV/c),l^^ or lower-px hadrons but only in minimum bias pPb 
events.l^^ Studies of the nuclear modification factor as a function of centrality in p- 
Pb collisions is highly challenging due to multiplicity selection biases, which remain 
inconclusive .1^^ Many studies have also indicated non-trivial correlations between 
the production of hard processes and underlying event activity, which is used for 
centrality determination in p-Pb collisions Qq clear path forward is 

laid out for the study of jet quenching in p-p and p-A collisions. 

Looking for azimuthal anisotropies (^ 2 ) of high-px particles can avoid the mul¬ 
tiplicity biases and provide us a hint of the (non-)existence of jet quenching. The 
fact that the near-side ridge yield persists to the high-px region suggests finite V 2 
values of high-px particles, as shown in Fig. 41 by the ATLAS collaboration.^^ 


Assuming back-to-back jet correlations are not significantly modified, the ATLAS 
collaboration extracted V 2 up to high px in high-multiplicity p-Pb collisions, where 
a sizable V 2 of about 5% is observed at px ^ lOGeV/c, shown in Fig. 41 

In Pb-Pb collisions, a long-range near-side correlation structure for a trigger 
particle with px > 20GeV/c is clearly visible (Fig. leftFurthermore, the V 2 
values obtained at high px and low px from different centrality ranges of Pb-Pb 
collisions are found to be strongly correlated (Fig. right), indicating a com¬ 
mon origin, i.e., that both are related to the initial-state geometry. It will be very 
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superSONIC: p,d,He +Au @ 200 GeV 



pt [GeV] 


Fig. 43. Left: The radial velocit y Vp extracted as a function of r (time) and p (radius) for a 
representative simulation from Ref.l^^with pre-equilibrium dynamics. Right: Ratios of ns between 
^He-Au, p-Au and d-Au as a functio n of p x at = 200 GeV from the superSONIC model 

with and without pre-equilibrium flow.E^ 


interesting to perform the same study in p-p and p-Pb collisions to gain some in¬ 
sight into the nature of long-range correlations at high px- Moreover, performing 
multi-particle V 2 measurements in the high-px region in p-p and p-Pb would help 
minimize the influence of short-range jet correlations and help clarify the picture 
of high-px azimuthal anisotropy. 


5.4. Pre-equilibrium dynamics 

The strong anisotropic flow of final-state particles in heavy-ion collisions has been 
regarded as evidence for thermalization realized during the very early stages of the 
collision. However, understanding the detailed mechanism of the fast thermalization 
process has been a big challenge in the field.In most hydrodynamic models, a 
QGP thermalization or formation time, tq, is usually assumed, after which the 
hydrodynamic evolution is switched on. This additional freedom allows for the 
tuning of models to match experimental data. Often, physics in the pre-equilibrium 
stage as well as its transition to a strongly-coupled QGP state is ignored, even 
though it is likely to have relevant effects on final-state observables (for example, 
its influence on HBT radii has been argued to be importanfP^. 

The IP-Glasma modeiMlIIIll includes some pre-equilibrium dynamics via the 
Yang Mills equations from time zero. However, a switch to hydrodynamics at a 
time To is still necessary, because the 2+1 dimensional Yang Mills simulation does 
not contain the necessary dynamics, which possibly includes the development of 
Weibel instabilities, to isotropize the system to the point where a matching with 
hydrodynamics can be completed. However, the build up of flow within the first 
^ 0.4fm/c is comparable in the Yang Mills and hydrodynamic picture, making the 
final state observables independent of the switching time tq.^ 

Progress has also been made in the superSONIG model, where pre-equilibrium 
dynamics and its transition to the hydrodynamic regime is treated in a coherent 
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way without any explicit! 


43 


(left) shows the development 


of the radial flow velocity profile as a function of time and radius going from the 
pre-equilibrium stage to the viscous hydrodynamic regimeP^ Most importantly, 
simulation results from the superSONIC model indicate that flow developed during 
the pre-equilibrium stage has a significant contribution for small colliding systems 
since the total lifetime of a small system is much shorter. This can be clearly seen 
in Fig. [ 4 ^ (right) for the ratios of U 3 between ^He-Au, p-Au and d-Au as a function 
of pt at = 200 GeV, compared with and without pre-equilibrium flow.^^^ 

A lot of work still lies ahead before any solid conclusion can be drawn but high- 
multiplicity, small colliding systems promise to open a new window for probing the 
pre-equilibrium dynamics in heavy-ion collisions. 


6. Summary 

In summary, the observation of long-range correlations in high-multiplicity p-p and 
p-Pb collisions has opened up new opportunities for investigating novel high-density 
QCD phenomena in small colliding systems. Experimental results from RHIC and 
the LHC over the past several years have provided crucial insights and imposed 
stringent constraints on possible theoretical interpretations. 

We reviewed the theoretical progress in the hydrodynamic modeling of small 
colliding systems. The current status is that hydrodynamic models are able to 
describe all the features of the experimental data at a largely quantitative level. 
However, a large sensitivity to the unknown initial state in small systems remains. 
Without further constraints on initial state models it will be difficult to distill infor¬ 
mation on the medium produced in small systems and scrutinize the hydrodynamic 
paradigm. Clearly hydrodynamics is being pushed to the edge of its validity in these 
small colliding systems and further theoretical work is needed to asses higher order 
corrections. 

Calculations based on the Color Glass Condensate effective theory can also de¬ 
scribe many characteristic features of the experimental data. If one can conclusively 
show that the correlations are generated by glasma graphs, valuable information 
about multi-gluon correlations in the nuclear wave-function can be obtained; multi¬ 
particle correlations can serve as a sensitive probe of saturation dynamics. 

Ideally these frameworks should be merged. We have a wealth of information 
from deep inelastic scattering that should be incorporated into our modeling of 
the initial state. How the initial gluon fields decohere, isotropize and possibly ther- 
malize to a system that can be described hydro dynamically is an open question. 
We have argued above that classical Yang-Mills dynamics contains many of the 
missing pieces of the glasma graph framework. Furthermore, the recently proposed 
color domain models can accommodate a negative C2{4}, and is highly interesting 
from a theoretical point of view as it could provide a measure of non-Gaussianities 
in the hadronic wavefunction. 

The strongest conclusion that can be made at this time is that the discovery 
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potential from small colliding systems is immense. The detailed experimental in¬ 
formation that has and will continue to come from these experiments will allow 
the community to test theoretical proposals at an unprecedented level of accuracy. 
Studying collectivity in high-energy proton-proton and proton-nucleus collisions 
provide us with access to a rich variety of emergent (and possibly yet undiscovered) 
QCD phenomena. 
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